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We suggest an approximative solution of the two dimensional Hubbard model close to half filling. It 
is based on partial bosonization, supplemented by an investigation of the functional renormalization 
group flow. The inclusion of the fermionic and bosonic fluctuations leads in lowest order to agreement 
with Hartree-Fock result or Schwinger-Dyson equation and cures the mean field ambiguity. We 
compute the temperature dependence of the antiferromagnetic order parameter and the gap below 
the critical temperature. We argue that the Mermin- Wagner theorem is not practically applicable 
for the spontaneous breaking of the continuous spin symmetry in the antiferromagnetic state of the 
Hubbard model. The long distance behavior close to and below the critical temperature is governed 
by the renormalisation flow for the effective interactions of composite Goldstone bosons and deviates 
strongly from the Hartree-Fock result. 
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I. INTRODUCTION 

The rich behaviour of strongly coupled electron sys- 
tems has attracted the attention of physicists for many 
decades. One of the most popular models, which is be- 
lieved to contain the most relevant features to explain as 
different phenomena as metal-insulator transitions, ferro- 
magnetism, antifcrromagnetism and, more recently, the 
behaviour of high temperature superconductors, is the 
Hubbard model Ij. It combines formal simplicity with 
a rich and complex spectrum of its implications that — 
even after nearly 40 years — are not fully understood. 
Analytical results are scarcely available. Only in one di- 
mension or in higher dimensions for very special choices 
of the parameters exact analytical solutions of the model 
are known Q . Besides many other attempts the last years 
witnessed strong efforts to investigate the properties of 
the model by means of renormalisation group equations 
0, 0, 0, 111 ■ Currently, although far from being able to 
fully and uniquely describe the properties of the Hubbard 
model, renormalisation group approaches are widely be- 
lieved to be a promising method to further investigate its 
interesting phase diagram. 

Most of the renormalisation group approaches have 
studied the evolution of the four fermion interaction in 
dependence on a renormalisation scale k. Phase transi- 
tions and spontaneous symmetry breaking are regarded 
to be encoded in the momentum structure of this ver- 
tex. More precisely, a phase transition is inferred from 
a divergence of the running coupling at some character- 
istic scale kssB- The particular momentum structure of 
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the divergent part of the vertex indicates the relevant or- 
der parameter for the phase transition. Even though the 
detection of a transition by this method looks rather con- 
vincing, the appearance of a divergence at /cssb forbids 
a continuation of the flow for k < /cssb- In consequence, 
it is not possible to penetrate the low temperature phase 
with non vanishing order parameter by a straightforward 
study of the four fermion interaction. 

For the two dimensional model this shortcoming is par- 
ticularly annoying since an exact theorem by Mermin 
and Wagner 7] forbids the spontaneous breaking of a 
continuous symmetry at non zero temperature. However, 
the widely believed antifcrromagnetism in the low tem- 
perature phase of the two dimensional Hubbard model 
corresponds precisely to the spontaneous breaking of the 
global spin rotation symmetry. Antiferromagnetic order 
at nonzero temperature therefore seems to contradict the 
exact theorem. The present paper will resolve this dis- 
crepancy and presents a computation of the antiferro- 
magnetic order parameter for temperatures below the 
critical temperature. A summary of our results can be 
found in Q. 

Our method is based on partial bosonisation, whereby 
the fermion bilinear ^ip is treated as a (composite) 
bosonic field (j>. For the Hubbard model, partial bosoni- 
sation replaces the four fermion interaction ~ U (ipip)'^ by 
a Yukawa coupling to the boson ^ h('ip'ip)(j) and a boson 

mass term ~ rn'^cjp, whereby U ^ h /ni^ . The renormal- 
isation flow describes how the mass term (k) depends 
on the renormalisation scale k. The onset of spontaneous 
symmetry breaking at fc = fcssB is signalled by the van- 
ishing of m^(fcssB)- This indeed corresponds to a diver- 
gent four fermion coupling U . However, in the partially 
bosonised formulation it is comparatively easy to include 
into the truncation also the bosonic interactions gener- 
ated by the flow, e.g. a quartic scalar interaction ~ At/)**. 
(In the fermionic language this corresponds to an eight- 
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fermion interaction which is quite difficult to compute.) 
For positive A a vanishing mass term = 0, or even 
a negative mass term < 0, still lead to a well be- 
haved effective scalar potential (or free energy density) 
UeS = rn^cjp /2 + A^^'/S. The minimum of the poten- 
tial shifts for rfi^ < simply to a nonvanishing value 
(pQ = — 2m^/A. The value of 0o can be associated with 
the order parameter responsible for spontaneous sym- 
metry breaking. For temperatures T below the critical 
temperature one finds that 4>{k) starts being different 
from zero precisely for k — ksgg > 0. It becomes obvious 
that the partially bosonised version of the renormalisa- 
tion group equation has no problem to describe the flow 
also for k < /cssb and therefore to penetrate the low tem- 
perature phase with spontaneous symmetry breaking. 

This allows us to study explicitely the issue of the 
Mermin- Wagner theorem. Our findings are indeed con- 
sistent with the formal correctness of this theorem. In 
an infinitely extended space all fluctuations arc included 
only if the effective infrared cutoff k is sent to zero. In- 
deed, we observe that (j)o{k ^ 0) = such that in the 
formal limit fc — > the order parameter vanishes and 
the Mermin- Wagner theorem is obeyed. However, in the 
low temperature phase below the critical temperature Tc 
the "running" of the k dependent order parameter to- 
wards zero is extremely slow - it is only logarithmic. As 
a consequence we find antiferromagnetic order at all real- 
istic macroscopic scales. We conclude that the Mermin- 
Wagner theorem is not practically applicable here, de- 
spite its formal correctness. 

We demonstrate this explicitely by computing the an- 
tiferromagnetic order parameter n and gap at a 
macroscopic scale — 1cm. The gap is directly con- 
nected to a nonvanishing antiferromagnetic order param- 
eter ao ~ <Pl by Aa = (27^^ao)^/^ = {2hlK)'^/'^t. Here k 
is the dimensionless renormalized order papameter ^ and 
ha the dimensionless renormalized Yukawa coupling such 
that the gap is proportional to the hopping parameter 
t of the Hubbard model. Fig^ shows the temperature 
dependence of the antiferromagnetic order parameter k 
at half filling (vanishing chemical potential n — 0). This 
figure may be regarded as a central result of our paper. 
We clearly observe the onset of spontaneous symmetry 
breaking for T <Tc- We find that the phase transition is 
continuous. For this plot we have concentrated only on 
the simplest version of the Hubbard model which includes 
only next neighbour hopping with coupling strength t. 
For a realistic value t = 0.3eV the critical temperature 
corresponds to Tc — 400K. 

An investigation of the antiferromagnetic phase within 
the partially bosonised formulation of the Hubbard 




^ Of course, in a finite volume V ^ the usual finite volume 
eff'ects also formally lead to a vanishing order parameter. This is 
common to all other systems with spontaneous symmetry break- 
ing. It is a well understood effect of no relevance to our issue. 



FIG. 1: Temperature dependence of the antiferromagnetic or- 
der parameter. We plot the renormalized expectation value k 
of the antiferromagnetic condensate in the low temperature 
phase for U/t — 3. The scale k/t = 10~^ corresponds to a 
probe of macroscopic size of roughly 1 cm. The gap for the 
fermionic excitations is related to n by Aa/t = \/27iJk. 



model has already been performed in Q where we used 
mean field theory (MFT). Unfortunately, we have found 
that in the MFT approximation the results for Tc or the 
antiferromagnetic order parameter depend strongly on 
unphysical parameters which describe the precise choice 
of the mean field. Since this ambiguity is closely related 
to the possibility of the Fierz reordering of the local 
four fermion interaction it is dubbed "Fierz-ambiguity" . 
In the present paper we show that the inclusion of the 
bosonic fluctuations in the renormalisation flow (beyond 
MFT!) greatly reduces this Fierz-ambiguity. 

In our present truncation the resolution of the mo- 
mentum dependence of the effective four fermion vertex 
depends on the choice of composite bosons included ex- 
plicitely. More precisely, the effective vertex arises from 
the exchange of the bosons and therefore inherits its mo- 
mentum structure from the momentum dependence of 
the bosonic propagators. In this paper we concentrate on 
the antiferromagnetic order and small doping. We only 
take into account bosons describing spin densities and 
charge densities. We emphasize, however, that our for- 
malism is quite general. It is also well adapted to include 
the interesting superconducting channel. 

This paper is organised as follows: In sect.^]we review 
the functional integral description of the Hubbard model 
and set up our notation. We proceed in sects. IIIII and 
IIVI to derive and solve the Schwinger-Dyson gap equa- 
tion. In sect. we introduce our partial bosonisation. 
We carry out a mean field calculation in sect. IVII that 
can be directly compared to the Schwinger-Dyson result. 
The "Fierz-ambiguity" arising in this approach is de- 
scribed in sect. IVIII In sect. IVIIII we turn to the mo- 
mentum dependence of the effective bosonic propagators 
which is first computed in the mean field approximation. 
The one loop corrections to the bosonic propagators show 
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that indeed antiferromagnetism is the dominant way of 
breaking the 5't/(2)-spin-symnietry of the model. The 
last part contains our renormalisation group approach 
to the Hubbard model in its partially bosonised version 
and particularly a discussion of the behaviour of the flow 
in the phase of broken symmetry. We review the exact 
functional renormalisation group equation for the effec- 
tive average action in sect. IIXI and the rebosonisation of 
fermionic interactions in sect. |3 Sect. IXII sets up our 
truncation and choice of cutoff functions for the Hub- 
bard model. In sect. IXlJ we derive the non-perturbative 
flow equations at half filling (/i = 0). Their solution in 
sect. IXHH leads to an understanding of the appearance 
of antiferromagnetic order despite the Mermin- Wagner 
theorem and is the basis for our main result shown in 
fig-D The critical behavior is discussed in more detail in 
section IXIVI In sect. IXVl we turn to nonzero doping and 
discuss the phase diagram of the Hubbard model for small 
values of the chemical potential. Section IXVII finallv dis- 
cusses how the Fierz ambiguity is substantially reduced 
if the effects of bosonic fluctuations are included beyond 
mean field theory. Our conclusions are summarised in 
sect. IXVIII In the appendices we give explicit formulae 
for the flow of the couplings in different truncations and 
show how nonlocal operators corresponding to a d-wave 
superconducting state may be incorporated into the for- 
malism. 



by 



II. THE HUBBARD MODEL 



The partition function of the Hubbard model is given 



m - / v{r{T)Mr)) 

J^(/3) = -^0(O),^-(/3) = -^-(O) 

+ \uY,Ui[T)Ur) 



(1) 



Here 



V'.W^d^jl^j), i^l{r) = {r.,{r\ r^.ir)) (2) 

are Grassmann fields describing electrons on a quadratic 
lattice with lattice sites at cc G Z x Z. (We restrict our at- 
tention to the two dimensional case.) The Euclideon time 
T parameterises a torus with circumference (3 = . The 
functional integral is constrained by antiperiodic bound- 
ary conditions due to the fermionic character of the fields. 



The kinetic term contains the chemical potential /i as well 
as the so called hopping matrix T, which describes the 
possibility of electron hopping between different lattice 
sites. For simplicity, we assume 



T — 



— t , if a; and y are nearest neighbors 
, else. 



(3) 



The quartic term reflects the local Coulomb repulsion of 
electrons at the same lattice site. Additionally we include 
fermionic sources rj^ir), rixir). 

Calculations are conveniently carried out in momen- 
tum space. We define the Fourier transforms 



Q 



rix) 



Q 



-iQX 



where we use the short hand notations. 

Q={uF,q), X = {t,x), QX = ujft 



xq 



(4) 



(5) 



in position and momentum space and define sums and 
delta functions accordingly 



E 



/3 



OG 



T 



cPq 



,^ ^-E' E-^ E J (2.) 

d{Q-Q')=(3Sr,^r.'{2TTf6{q-q') 

5{X-X') = 5{t-t')5^,^,. 



(6) 



Here uip = {2n + 1)-kT is called Matsubara frequency and 
the corresponding sum Matsubara sum, with inverse tem- 
perature f3 = T^^. All components of Q or A" aquire the 
same canonical dimension if q and x are measured in 
units of the lattice distance, i.e. x ax, q — > q/a. Here 
we set a — I. Note that all functions in momentum space 
are 27r-periodic in spatial momentum. In particular, the 
kinetic term in momentum space takes the form 



S 



F,ki„ = E^^(Q)^(Q)^(Q)' 
Q 



(7) 



P{Q) — iujp ~ — 2t(cosqi + cos 52), 



such that the Fermi surface is located at —2t{cosqi + 
cos 52) = ^J■■ 

Another convenient notation that we will adopt is to 
introduce a collective index a denoting spin, A or Q, and 
the label for the two complex conjugated fermionic fields. 
Correspondingly, a = 1 ... 4 is the collective index associ- 
ated with the discrete "internal" indices only. Explicitly, 
we have 



MX) 



(X). 



In this notation, the partition function reads 
Z[ri] = / Vipa exp (^-S[ip] + rja-ii^o^ 



(8) 



(9) 
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with the action 

= -^ia/3?/'a# + ^KpiS^ai'pi^^^S- (10) 

In momentum space jap — iab{P,Q) is defined by the 
Fourier transforms implying 



jab{P,Q)=JadQ)5{P-Q). 

P(-Q)l2 



J{Q) 



-P{Q)i2 



(11) 



The quartic coupling Xap-fS is most easily given in posi- 
tion space. By comparing ij^ with we get 

Kp^s - -t/eabcd 6{X - X2)S{X - X^)5{X - X^) (12) 

with a = {X,a), (3 — {X2,b), 7 — (^3,0) and S = 

{X4,d). 



III. SCHWINGER-DYSON EQUATION 

Starting from the partition function we derive 
the Schwinger-Dyson equation [lol | for the fermionic two 
point function. The lowest order approximation to the 
Sch wing er-Dyson equation is the Hartree-Fock equation 
[Til IT a IT^ [3 . We start by defining the effective ac- 
tion by the Legendre transform of the logarithm of the 
partition function in the standard way: 



r[V'] = -lnZ[77[^]]+ry„[V^]V'a. 



(13) 



Here -tpa = {''Pa) — denotes the "expectation 

value" of the corresponding fermionic field in presence of 
fermionic sources. It is possible to rewrite the effective 
action as an implicit functional integral 



r[i^] = - exp + ^] + 77,^;,] (14) 

with rja — —dT/dipa- We will use this as the starting 
point for the derivation of the Schwinger-Dyson equation 
[l5'|. For vanishing sources one finds the derivative of eq. 

m 

dThp] 1 ... 

= -jaai^a + ^Ka/Sy ii'a'^ P^Jj) (15) 

+(r(^^);>7-(r(^^);^.V^ + (r(^^),^,V. 



where 



dipadipp 



(16) 



is the inverse fermionic propagator. Equation (|15|) can 
be read as an identity relating various fermionic n-point 



functions. Differentiating (|15|l once more and setting ip 
0, we find the Schwinger-Dyson equation 



^ per ^ Jpcr ^ cy'^pcraPK^ )af3 



(17) 



relating the fermionic 2- and 4-point functions. The last 
term in this equation represents a two loop expression 
of the order j3| and will be neglected in this and the 
next section. The lowest order Schwinger-Dyson equation 
results in a closed equation 



1 



(18) 



— —ja(3 + -^^apjsGjS 

for the two point function Gap, 

Gap = (r(^));^. (19) 

It can be integrated in order to obtain the two particle 
irreducible effective action or bosonic effective action , 
|15| and a corresponding expression for the free energy. 

IV. ANTIFERROMAGNETIC GAP EQUATION 

The Hubbard model is symmetric under global 
5?7(2)-spin - and [/(l)-charge-transformations. However, 
for small T and small /i (corresponding to small doping) 
it is well known that the Hubbard model exhibits antifer- 
romagnetic behaviour. A nonvanishing expectation value 
of the operator ip'^ {X)(jip{X) in macroscopically large re- 
gions can be associated with spontaneous breaking of the 
S'C/(2)-spin-symmetry. Even though spontaneous sym- 
metry breaking of a continuous symmetry is forbidden in 
a strict sense by the Mermin- Wagner theorem Q we will 
see in section lXIIII that the main physical features are cor- 
rectly described by such a "naive spontaneous symmetry 
breaking" . The leading order Schwinger-Dyson equation 
exhibits a nonzero expectation value^ of the antiferro- 
magnetic condensate^ 



9ab 



U 



■ exp(iHX) 



-a 
?^ 



(20) 



with n = (0, TT, tt). Here the factor exp(iHX) ensures the 
alternating sign of the spin density between neighbouring 
lattice sites. We therefore investigate the possible solu- 
tions of equation 1)18(1 for Green functions obeying 



Gab[X,X)^g'^^{X) + 



(21) 



^ Note that the result of the leading order Schwinger-Dyson equa- 
tion is not consistent with the Mermin-Wagner theorem. Our 
more extended approximations based on the renormalisation 
group equation will establish the missing consistency. 

^ The Yukawa coupling ha is introduced solely for consistency with 
later notation and plays no role here, only haS being observable. 
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where the symmetry conserving part can be parame- 
terised as 



hpP f -I2 
U \l2 



(22) 



Since the second term on the r.h.s. of equation H18|) is 
local (due to equation dJ) it only involves GabiX,X) 
and therefore induces a "local gap" 

Aafc(X) = -^^abM^\X)+g^^\X)) ^^^^ 
= (Ap)afc + (Aa)ab{X) 

with charge density shift 

'0 -I2 



Ap = -hpp 



I2 



and antiferromagnetic spin density gap 

(Ka)ab{X) = {Aa)abexp{mX), 



Aa = -ha 



-a^ a 
ad 



(24) 



(25) 



For nonvanishing Aa the fermionic propagator indeed ac- 
quires a mass gap according to 

G-,\X,Y) - -jabiX,Y) + AabiX)S{X - Y). (26) 

On the other hand, we will see below that the charge 
density shift only shifts the effective value of the chemical 
potential, Me// = M = ^pPif^) with p(/i = 0) = 0. 

In order to determine had and hpp as functions of T 
and fi we have to solve the "gap equation" which ob- 
tains by evaluating the inverse of equation H26() at equal 
arguments X = Y, i.e. 

g{X) = GiX,X)^{-j+A)-\X,X). (27) 

In momentum space this becomes 

g{Q) = J2 G{K, Q + K) = ^(-j + A)-\K, Q + K) 

(28) 



K 



K 



with gap function 

A{P,Q)^Ap5{P~Q) + Aa6{P-Q + n). (29) 



The solution of (|28|l needs the inversion of the inverse 
propagator H26() in momentum space 

(-j+A)(P,Q) = {-J{P)+Ap)S{P~Q)+Aa&{P-Q+'R)- 

(30) 

By use of the identity 

A:(-J(Q) + A,)* = (-J(Q) + Ap)Aa (31) 

we find 

(-j+A)-i(P,Q) = 

= -N-\P){{J{P + Ii)-ApY5{P-Q) (32) 
+ Al5{P + n-Q)] 



with 

N{P) = iJ{P + n) - Ap)*(J(P) - A,) - AlAa 

= A(P) + 2it^FAicff diag(l, 1, -1, -1), 
A{P) :^uj% + hld^ - nlff + 4t2(cospi -I- cosp2)^ 

Inserting eq. H32() into eq. (|28l) and performing the Fourier 
transform of eq. ^ (with g{X) J2q e''^^ 9iQ)) we 
arrive at the gap equation for had 



A{Q) 



(34) 



Equation H34|) admits a solution with a ^ for low 
enough T and /Xeff. This nontrivial solution has indeed a 
lower value of the relevant free energy ^3 compared to 
the solution with a = 0. It is therefore a good candidate 
for the thermal equilibrium state at low T. (We will argue 
in more detail below (sect. VIII) that this state indeed 
minimizes the appropriate free energy in the approxima- 
tion of the lowest order Schwinger-Dyson equation.) The 
onset of spontaneous symmetry breaking (a ^ 0) ~ as the 
temperature is lowered from high values where a = - 
can be associated with a phase transition. We have solved 
eq. (|34ll numerically P| and the result for the critical line 
in the plane (T, pcs) is shown in figI21 (upper curve) . 




FIG. 2: Phase boundary for low temperature antiferromag- 
netism in the lowest order Schwinger-Dyson and mean field 
approximations. We show the critical temperature as a func- 
tion of the chemical potential fi. We choose couplings {Um = 
(7/3, Up = 0) (lower line) and {Um = U/2, Up = (7/2) (up- 
per line). The upper line corresponds to the the solution of 
the Schwinger-Dyson equation. Both curves correspond to the 
same value U/t — 3 and their difference reflects the Fierz am- 
biguity discussed in sect. IVlfl 

We can extend the analysis for other possible forms 
of the local gap (or g{X)). For example, ferromagnetism 
(instead of antiferromagnetism) can be described by re- 
placing in eq. H2U|) hadexpiillX) hmrh. We have found 
that the solutions with m 7^ have a higher free en- 
ergy as compared to a 7^ and that antiferromagnetism 
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is indeed favoured. Charge density waves (with an ad- 
ditional factor exp(inX) in eq. H22I) ') do not aUow for 
nontrivial solutions. The same is true for local Cooper 
pairs with s— , p— or d— symmetry. Nonlocal symmetry 
breaking expectation values Gab{X, Y ^ X) do not con- 
tribute to the lowest order Schwinger-Dyson equation 
due to the locality of the four fermion interaction. Within 
the Schwinger-Dyson approach an investigation of real- 
istic nonlocal Cooper pairs which could be responsible 
for superconductivity therefore requires the inclusion of 
the higher order terms. (See [l5 ] for an explicit formula 
of the next order term.) As an alternative, one may re- 
alize that the renormalisation flow of the four fermion 
interaction introduces nonlocalities in the effective four 
fermion interaction Q and use the "renormalized gap 
equation" proposed in |l5j| . For nonlocal four fermion in- 
teractions even the lowest order Schwinger-Dyson equa- 
tion can have nontrivial solutions for nonlocal gaps. 

We will go beyond the leading order Schwinger-Dyson 
equation by using renormalisation group equations from 
section lDO on. Before, it will be instructive to understand 
the connection between the Schwinger-Dyson equation 
and mean field theory based on partial bosonisation. This 
will provide us with an appropriate setting for the renor- 
malisation group equation by which we can explore the 
region of spontaneous symmetry breaking. 



V. PARTIAL BOSONISATION OF THE 
HUBBARD MODEL 

Among the most promising approaches to investigate 
the properties of the Hubbard model are renormalisation 
group studies. So far, these calculations have been mainly 
carried out in a purely fermionic formulation U, |5j • 
The onset of spontaneous symmetry breaking then shows 
up in the divergence of the quartic coupling in certain 
momentum channels, if one follows the renormalisation 
group flow to larger length scales. This approach has 
the limitation that one is not able to follow the flow 
into the region of broken symmetry. This would need the 
inclusion of higher multifermion interactions and there- 
fore become extremely complex. We adopt here an alter- 
native approach, where symmetry breaking is encoded 
in non vanishing expectation values of bosonic fields. 
Each bosonic field corresponds to one possible symme- 
try breaking channel. This allows us to keep a simple 
overview in the realistic case of competing channels. In 
this work we will work with a bosonic field reflecting an- 
tiferromagnetic behaviour. Similarly, bosons correspond- 
ing to S-, p- or d-wave superconducting behaviour could 
be defined for an extended discussion '9'|. 

The main idea is to rewrite the original Hubbard model 
with a quartic fermionic interaction as an equivalent 
Yukawa theory, where the quartic fermionic interaction is 
replaced by fermion bilinears coupled to the bosons [l8j |. 
As mentioned above, symmetry breaking shows up as 
non vanishing expectation values of some bosonic fields. 



These order parameters can be computed by following 
the renormalisation group flow for the bosonic effective 
potential (free energy density) into the broken symme- 
try region. In particular, it is now straightforward to 
investigate the role of bosonic self interactions which 
correspond to higher multifermion interactions in the 
purely fermionic language. Furthermore, the identifica- 
tion of symmetry breaking channels with bosonic excita- 
tions (particles) enhances our intuitive understanding of 
what happens during the flow. 

However, this approach, as advantageous it may seem 
at flrst sight, is plagued by a coupling ambiguity prob- 
lem. In the fermionic theory, we face one single quartic 
coupling U, whereas in the partially bosonised theory, 
we introduce as many Yukawa couplings as we introduce 
bosonic fields. These Yukawa couplings are not uniquely 
determined by the equivalence to the original fermionic 
model - different choices of these couplings lead to pre- 
cisely the same fermionic model if we integrate out the 
bosons. This coupling ambiguity would not give rise to 
any problems were we able to solve the flow equations 
exactly. But for many approximations the results will un- 
physically depend on the choice of the couplings. At this 
point, the Schwinger-Dyson approach presented in the 
last sections comes in handy. We will see that at least 
for a given channel the gap equation corresponds to the 
mean field equation with a definite choice of the Yukawa 
couplings. This determines at least the rough range for 
the Yukawa couplings which can serve as a reasonable 
starting point for approximations beyond the mean field 
approximation. We will discuss in sect. XV how the in- 
clusion of the bosonic fluctuation effects by the use of the 
partially bosonized flow equations reduces drastically the 
dependence of physical quantities on the precise choice of 
the Yukawa couplings. 

In this section we will show how to partially bosonise 
the Hubbard model for the case where we only take into 
account a charge density boson p{X) and an spin density 
boson rn{X). These are the same degrees of freedom we 
used in the ansatz (|24() and H25() for the gap function. Our 
bosonisation procedure preserves explicitely the continu- 
ous S'?7(2)-spin-symmetry while it is consistent with the 
Hartree-Fock result. This settles a widely discussed 
apparent conflict betwen these two aims. We begin by 
defining the fermion bilinears 

p{x) = ^Hx)4'ix), 

, (35) 
m(A) = i;\X)a'tlj{X). 

With these bilinears the fermionic coupling term of the 
original Hubbard model can be written as 

(^t(x)V'(A))2 = piXf = -^rfiiXf. (36) 

Suppressing the explicit notation for X the fermionic ac- 
tion in the presence of sources for the bilinears is given 

by 

Sf = SpMn + ^Uiip^ip)"^ - Jpp - J„im, (37) 
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with SpMn the fermion kinetic term (j?)). The partition 
function reads 



Sn 



(38) 



We define the partially bosonised partition function by 



Z[7-i,ri*,Jp,Jjn]= J V{i{j*,'ilj,p,m)exp{-{S + Sr^ + Sj)) 

(39) 

with 

S = 5'F,km + ^Upp'^ + \Umrn^ ~ Uppp - Um-mm, 
Sj = -JpP - Jmm. 



(40) 



The functional integral now also involves the bosonic 
fields p and m. The equations of motion for the bosons 
yield 



P 



Uo 



m = m -\- 



(41) 



We now demand that inserting the bosonic equations of 
motion, S should reduce to Sp up to field independent 
terms. It is then easy to check by Gaussian integration 
that the partition function (|39|l is equal to H38|l up to 
the logarithm of a field independent quadratic function 
of the sources . Equivalence between the bosonic and 
fermionic formulation is therefore achieved for 



u = -Up + m„ 



(42) 



At this point we explicitly see that different combinations 
of the couplings Up and Um in the partially bosonised 
theory lead to the same fermionic theory. 

In summary, the fermion-boson action takes the fol- 
lowing form in momentum space 



S = S'f.Mii + Sb + Sy + Sj, 



(43) 



'F.kin 



' ip^Q){iuJF - M - 2t(cosgi + cosq2))tp(Q), 



= J 5] (t^pP(Q)/5(-Q) + M(Q)'^(-Q) 



Q 

Sy = ~ 6iQ-Q' + Q") X 
QQ'Q" 

{Upp{Q)i^\Q'MQ") + u,^va{Q)^\Q')'^i^{Q"\ 

Q 

Here the bosonic momenta Q — {ujB,qi,q2) involve Mat- 
subara frequencies ub — 27rmT with integer m. 



fields to constant background fields. Only the fermion 
fluctuations are included in the functional integral 
In this "mean field theory" (MFT) the fermionic integra- 
tion can be carried out analytically. In particular, con- 
stant bosonic "background fields" p and a 



(44) 



p{Q) ^ p5{Q) 

rh{Q) aS{Q~n), 

correspond to homogeneous charge density (p) and an- 
tiferromagnetism (a). The mean field partition function 
reads 



Zmf = J V{r,t^)eM-SMF), (45) 

Smf = ip'^{Q){iujF - p - 2t(cosgi + cosg2))'0(Q) 
Q 

- Y.iUpPi'HQMQ) + UmS^^iQ + n)a^(g)) 

Q 

+ ^iUpP^ + Urna^) - Jp{Q)p - J™(-n)a, 

with V2 the two dimensional volume of the lattice. From 
the effective action 

TmF = - In^MF + Jp(0)p + Jra{-^)d (46) 

we can derive the effective potentiate 
TV 

Uip, d) = — = ^{Upp^ + U„,a^) + AU{p, a), (47) 
y2 

with the fermionic one loop contribution 
T 



AUip,d) ^-y^^^J 2?(V'*,V')exp(-5A), (48) 
- ^ (v>(Q)P(g)V^(Q) - U„,ai;HQ + n)a7^(Q)) , 



and P{Q) now given by 

P{Q) — iujp - pcff - 2t(cosgi -I- cos 92), 
PcS = P + Upp. 
By defining 

we can cast Af/ into the form 
AC/(p,a) = 

= -^in j vir,i')exp |^-^^v^^(Q)i^(g)v;(Q)j 



(49) 



(50) 



P{Q) -U^dd 
-U^dd P{Q + U', 



(51) 



VI. A MEAN FIELD CALCULATION 



In this section we will explicitly calculate the effec- 
tive action in an approximation that sets the bosonic 



* There should be no confusion between the effective potential and 
the coupling in the Hubbard model which traditionally arc both 
denoted by the letter U. 
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The integral is Gaussian and can be carried out 

AU{p,a) = -i^ln(-detP(Q)) (52) 
Q 

= - ^In (^^1 - /i^fj + 4t^(cosqi + cos(j2)^ 
Q 

The Matsubara sum in AU can be analytically per- 
formed, so that we finally find (apart from a temperature 
dependent constant) 



AU 



(53) 



££{-1,1} 



2^ (^cff + eV4t2(cos(7i + cosgz)^ + U'^a^^ ^ . 



In order to establish the connection between the mean 
field approach and the lowest order Schwinger-Dyson 
equation we take a closer look at eq. (|52|l . The value 
of minimising the effective potential is zero for the 
symmetric phase, but differs from zero in the case that 
the the symmetry is spontaneously broken, which in our 
case means that the system shows antiferromagnetic be- 
haviour. The necessary condition for this is = 0, 
which can be brought into the form 



1 = 2Ur 



A{K) 



(54) 



This is the same as eq. if we identify the gap hpp 

in the gap equation Csect. lIV|) with Upp here (so that pcS 
used in the gap equation and the mean field calculation 
coincide) and similarly C/,^a^ = h^a^ ■ We also have to 
set Um = UI2, which imphes Up = U/2 by eq. lO- The 
Schwinger-Dyson approach therefore does not have the 
coupling ambiguity introduced in the partial bosonisa- 
tion procedure. As mentioned before, it corresponds to a 
special choice of the couplings in mean field theory. 



VII. MEAN FIELD AMBIGUITY 

For high enough temperature and small enough cou- 
pling an expansion in the small dimensionless quantity 
U /T is expected to be valid. In this limit the lowest or- 
der Schwinger-Dyson equation becomes the lowest order 
in a systematic expansion where higher orders are sup- 
pressed by higher orders of U /T . (This is not a simple 
expansion in powers of U /T like perturbation theory but 
rather corresponds to a resummed expansion.) For small 
U /T mean field theory therefore gives an incorrect re- 
sult unless we choose the particular bosonisation proce- 
dure with Um = Up = U/2. Since the functional integral 
transformations leading to partial bosonisation are exact 
the failure of mean field theory for other choices of Um 



and Up must be due to the bosonic fluctuations which 
have not yet been included. If we were able to compute 
the effects of the bosonic fluctuations exactly the flnal 
results would have to be completely independent of the 
choice of [/,„ provided eq. H42|l is obeyed. In general, the 
inclusion of the bosonic fluctuations is therefore crucial 
for a reliable picture. This issue has been systematically 
studied in 19]. For small U/T it just happens that the 
sign of the effect of the bosonic fluctuations may be pos- 
itive or negative, depending on the value of Um- In par- 
ticular, we infer that to lowest order in U/T the bosonic 
correction to the (anti)ferromagnetic correlation function 
{m{X)m{Y)) precisely vanishes for Um = U/2. 

Near the phase transition an expansion in small U is no 
longer reliable, the characteristic quantity U/\T—Tc\ even 
diverges for T —f Tc. The corrections to the lowest order 
Schwinger-Dyson equation become important. In conse- 
quence, there is no reason to believe that the bosonic 
fluctuations beyond mean fleld theory are a small effect 
for Um = U/2. We conclude that for the range of large 
U/\T — Tc\ the "gap equation" value Um = U/2 is at best 
a good educated guess. We feel that it is much safer to 
study a whole range^ of choices for Um- The dependence 
on the choice of Um/U in mean field theory may then 
be turned into an advantage when fluctuations beyond 
mean field theory are included: if an approximation be- 
yond mean field theory is valid the dependence on the 
unphysical parameter Um/U should get considerably re- 
duced as compared to mean field theory. This becomes a 
demanding check for the validity of approximations. We 
will see in section IXVII that our renormalisation group 
approach indeed has this important property. For the mo- 
ment, we illustrate the mean field ambiguity by the lower 
curve in the phase diagram in fig. [21 which corresponds 
to the choice Um = U/3, Up = 0. In the language of the 
present section this figure shows the mean field results 
for the onset of spontaneous symmetry breaking towards 
the antiferromagnetic state for two choices of the Yukawa 
couphngs Up,Um, namely [Um = U/3, Up = 0) (lower 
line) and for [Um = U/2, Up = U/2) (upper line). The 
lines indicate a second order phase transition and stop 
where the phase transition becomes of first order. We 
observe that without the inclusion of the bosonic fluctu- 
ations the "mean fleld ambiguity" for the phase diagram 
is considerable. Depending on the choice of Um/U the 
critical temperature varies by more than a factor of two, 
limiting severely the reliability of mean field theory. 



^ One should, however, avoid the extreme scenarios where Um/U 
becomes very large or very small or negative. A reasonable range 
may be within 30% of Um/U = 0.5. 
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VIII. MEAN FIELD CALCULATION OF THE 
BOSONIC PROPAGATOR 

In this section we will calculate the bosonic propagator 
for m in mean field theory. This also gives the bosonic 
propagator according to the lowest order Schwinger- 
Dyson equation if we choose Um — U/2. There are two 
reasons why we are interested in the result. In the deriva- 
tion of the gap equation and the mean field approxima- 
tion, we have assumed that antiferromagnetism is the 
dominating mechanism by which the SU (2)-symmetry 
is spontaneously broken. By inspection of the bosonic 
propagator as a function of momentum we will see that 
indeed antiferromagnetism (in contrast to e.g. ferromag- 
netism) is the favoured mode of symmetry breaking. The 
second reason is that — as we will see in the subsequent 
sections — our renormalisation group equations take the 
form of one loop equations. The calculation in this sec- 
tion therefore can be immediately used as a technical in- 
gredient for the derivation of the renormalisation group 
equations. Furthermore, the general features of the mo- 
mentum dependence of the bosonic propagators found in 
this section will be useful in order to device truncations 
for the flow equations. 

The bosonic propagators are propagators for fermion 
bilinears and therefore correspond to particular combina- 
tions of four fermion correlation functions in the original 
fermionic description. They are particularly simple to de- 
scribe in the partially bosonised language of section W\ 
For example, the propagator for the spin wave boson 

K(X)m,(y)), = (GS)»,(X,y) (55) 

is given by the inverse of the second functional derivative 
of the effective action with respect to the scalar fields in 
the appropriate channel. Here the effective action in the 
partially bosonised language obtains from the partition 
function by the usual Legendre transform 

r[p, m, i/i, -0*] = -lnZ + Jpp + JmTO + 77V'' + '7^0'* • (56) 

As usual the second functional derivative of F with re- 
spect to p and m is directly related to the second func- 
tional derivative of In Z with respect to the sources Jp 
and J„j (by inversion). Applying this to eqs. I|38|) and 
(|37|l and setting Jp — Jm = we see indeed how partic- 
ular combinations of fermionic four point functions can 
directly be read off from the bosonic inverse propagator 

m 

In absence of spontaneous spin-symmetry breaking the 
inverse bosonic propagator for to obtains by expanding F 
in second order in rfi for fixed constant p and ip — ip* = Q 

T^\Y,^^i~Q)i'^t%iQ,Q>lAQ') + ■■■ ^ (57) 

QQ' 

vZ\Q,Q') = V^^\Q)6{Q~Q') (58) 

(with G^'' ~ (Fm'')^^). In section VI we have already 
computed t'^m {Q — H) in equations lf3^ . It corre- 



sponds to the antiferromagnetic susceptibility or inverse 
correlation length 

(f (^)),,(n) = rn%„ ml = U^ + 2^AUy=o- 

. . " . 

The critical line for the phase transition shown in fig. 
121 corresponds to rn^{T = Tc) = 0. If we reexpress this 
in the fermionic language the four fermion interaction 
in the antiferromagnetic channel diverges for T — » Tc 
proportional to m^^. 

In this section we want to compute Fm (Q) for arbi- 
trary Q. Antiferromagnetism is the preferred condensate 

" (2) 

(with lowest free energ y) if FJ^'(Q) has its minimum for 
Q = n. (A minimum at Q = would favour a ferro- 
magnetic condensate.) We start from the action H43|) and 
treat the bosonic fields as background fields according to 

p{Q)-^P, rfi{Q)^m{Q). (60) 

Note that we now keep the momentum dependence of to 
to be able to distinguish between ferromagnetic, antifer- 
romagnetic, etc. behaviour. For a (mean field) computa- 
tion of the fermionic fluctuation effects we need the ex- 
pansion of the action 1431) in second order in the fermion 
fields 

S2^Y. ^^(Q)iPo{Q,Q') + AP(g,Q'))V^(Q') (6i) 

QQ' 

with 

Po{Q,Q') = P{Q)S{Q-Q') 
AP{Q,Q')^^Urnm{Q~Q')a ^ ' 

and P{Q) given by eq. (|43l) in second order in the fermion 
fields The fermionic functional integral can be carried out 
and we get 

AF = - In j I?(V3*,V')exp(-52) (63) 

= -lndct(Po + AP) 

= AFo + AFi + AF2 + • • • 

= -TrlnPo - Tr(Po"^AP) + ]^T:t{P^^ APP^^ AP). 

The first term is independent of m, the second vanishes, 
and AF2 is the TO-propagator correction we are looking 
for. Its second derivative is given by 

(Af(^)),,(Q) = AP„,(Q)% (64) 

AP™(g) = 2t/^^p-i(i^)p-i(i^ + g). 

K 

If we set g = n, we reproduce the term in the mean field 
result of quadratic order in a. 

The inverse propagator for the spin wave bosons can 
now be computed numerically after adding the classi- 
cal piece (gni) S'o, i.e. Pm(wB, 91,92) = fm^g) = Um + 
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/S.Pm{iOB, 9i, 92)- For vanishing Matsubara frequency (i.e. 
Q = (0, 52)) and for ^ = the Matsubara sums in 
can be performed analyticaUy 



AP„,(Q) = U„ 



2T r (fq 



Sim, 



Eq+fc - 



•(7rT)2 (27r)2"'^--'' ttT' ttT 
For this purpose we use the identity 

ab- {2n + l){2{n + m) + 1) 



(65) 



S{m;a,b) = ^ 



a2 + {2n + 1)2] [62 + (2(n + m) + 1)2] 

(66) 



TT (a - 6) (tanh ^ - tanh ^) 



2 ' 



(67) 



2 4m2 + (a - 6)2 

which has the following limits 

7r2 

5(0; a, a) = —— cosh " 

TT tanh if 
2 a 

We show the momentum dependence of P,„ in fig.Oand 
observe a pronounced minimum for Q — II. This clearly 
establishes that antiferromagnetism is the preferred con- 
densate as compared to other spin waves with Q ^ H. We 
also note that the minimum of P„i gets more pronounced 
as the temperature is lowered, reaching Pm(n.) = for 
T = Tc. (For our example Tc « 0.2i.) Also notice the 
development of sharp crests at low temperature (lower 
figure) due to the singularities in the fermionic propaga- 
tors at the Fermi surface. 

Fig.^shows the dependence of Pm{Q) on the Matsub- 
ara frequency ujb = 2'KmT for two values of the external 
momenta. Note that the = mode is the one that 
is changed most by fluctuation effects. For = q2 = Q 
this is even the only mode that is changed at all, while 
APm vanishes in the other cases. We observe that Pm is 
even in as requested by the discrete (time reversal) 
symmetry of the Hubbard model. 

A main reason for calculating the one loop corrections 
to the bosonic propagators was to get a feeling for the 
momentum dependence the propagators are likely to ob- 
tain under a renormalisation group flow. In sect. XI we 
will make use of these observations in order to formu- 
late suitable truncations for the bosonic propagators. A 
more detailed discussion of the propagators of the other 
bosonic fields (charge density waves, Cooper pairs) can 
be found in p7| . 



IX. EXACT RENORMALISATION GROUP 
EQUATION 

In the remainder of this paper we want to extend our 
analysis of the Hubbard model beyond the approxima- 
tion of mean field theory. Our tool will be nonperturba- 
tive flow equations for the scale dependence of the ef- 
fective average action [23. The flow equations are ob- 
tained from a truncation of the exact functional renor- 
malisation group equation (ERGE) |^ . Within the 





FIG. 3: Mean field approximation for the inverse propagator 
of spin waves as a function of momenta. We show two different 
temperatures T = 0.5t (upper part) and T = 0.15f (lower 
part) for /icff = 0, ujb ~ and choose Um ~ t. 



fermionic language a truncation to momentum de pen dent 
four fermion interactions for the ERGE of ref. [23| has 
been performed^ in 0. This analysis shows for T ^ Tc 
and /i = a leading divergence of the four fermion in- 
teraction in the antiferromagnetic channel. It also covers 
the issue of superconductivity for large enough and 
the analysis has been extended beyond next neighbour 
coupling. 

The limitation of the fermionic language concerns the 
region when some of the couplings grow very large. One 
expects that large quartic fermion couplings also induce 
large eight fermion couplings and so on. In particular, 
it seems very difficult to explore the low temperature 
phase with spontaneous symmetry breaking in this way. 
For this reason we explore here the formulation in the 
partially bosonised language of section As we have 



® Ref. 1^ uses a somewhat different formulation of the exact renor- 
malisation group equation. 
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rv^ 0.7 



0.7 - 
0.65 - 



FIG. 4: Inverse propagator for the boson m as a function of 
the Matsubara frequency ujb = 2nmT a± T = 0.2t for Um = t- 
We show gi = 92 = (upper part) and qi = q2 ~ tt (lower 
part). 



surface (fermionic fields). In presence of the cutoff the 
generating functional for the connected Green functions 
depends on the scale k. 



Wk[J] 



In I Vxexp{^{S[x]+ASk[x]) 



JX)). (70) 



More precisely, the function Rf{Q) is tailored in order 
to regularise the fermionic zero modes of the propagator. 
For momenta far from the Fermi surface (compared to 
k) Rf is to vanish rapidly so that the behaviour of these 
modes is essentially unaltered. A similar task is assigned 
to the bosonic cutoff functions. Here the zero modes can 
occur for = on the critical line and i?^!(Q), K^{Q) 
act essentially as additional mass terms ^ k^ for the low 
momentum region. This is easily generalised to the case 
when the zero mode occurs for some other momentum, 
e.g. Q = n. In the limit A: ^ we demand that the 
regulators vanish whereas for A: — > cx) we assume them to 
diverge 



lim = 

fc^O 



0, 



lim R 



(71) 



Our specific choice of the cutoff functions will be given 
in sect. IXII 

We may now proceed to define the effective average 
action F^, in analogy to the definition H13|) . By a Legendre 
transform with respect to the classical fields 



we obtain the functional 



SJ 



(72) 



argued, the divergence of a four fermion coupling just 
corresponds to the vanishing of the mass of the bosonic 
field. If higher order bosonic interactions are included 
in the effective bosonic potential U one can easily ex- 
plore the region where the curvature of U at the origin 
gets negative (m^ < 0). This region is characterised by 
"spontaneous symmetry breaking" . 

Let us consider a theory containing complex bosonic 
fields u, u* , real bosonic fields w and fermionic fields ip, 
ip*. We collect the fields into generalised fields and define 
generalised sources for them (the indices a run over field 
type, momentum, internal indices etc.) 

Xa = {u,U*,W,1p,1p*)a, (68) 
Sj[x] = -JaXa = -{.fu + ju* +lw + ■q*'4) + 771/^*). 

Now we regularise the theory by adding an infrared 
cutoff ^ Rk to the original action 

ASk[x] = kxcRKcpXf) - u*Rlu + \w*R^w + rRU- 

(69) 

It will have the effect of cutting off momenta with < fc^ 
(bosonic fields) or in a shell of size k around the Fermi 



fk[x]^Jx~Wk[Ji 



(73) 



where J — J[x] is a solution of the equation H72|l . As will 
become clear in a moment it is favourable to subtract the 
cutoff action from this functional and define the effective 
average action as 



Tk[x]^Jx^Wk[J]-ASk[x]- 



(74) 



The effective average action is the effective action of a 
theory containing an extra regulator or "mass" term as 
described by the action S[x] + ASk [x] ■ Since the effective 
action respects all (linearly realized symmetries of the 
original action, this also applies to Tk[x] for all k, if the 
regulator AS'fe[x] respects the symmetries. It is thus pos- 
sible to expand the effective average action in invariants 
with respect to these symmetries. 

The limits (|71|l lead to corresponding limits^ for the 
effective average action |23| 



limr.M 



nx], 



lim Tk[x] ^ S[x]. 



(75) 



For Rk 
stant). 



00 the fluctuation contribution vanisiies (up to a con- 
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This is why we chose to subtract the regulator in the 
definition of rfc[x]: for large "cutoff" k this functional 
is nothing but the original action. The effective average 
action therefore interpolates smoothly between the clas- 
sical action S and the effective action F as fc changes from 
very large values to zero. The quantitative change of Tk 
with k will be described by the ERGE derived below. 

We specify the second functional derivative in a sym- 
metric form containing both left and right derivatives 

and notice the identity 

^klMX = M^^, M = diag(l, 1, 1, -1, -1). 

(77) 

For the fc-derivative of f ^ one now obtains 0, : 

dkTkMx = -dkWk[J]\j = {dkASk) 

^^d,Rk,o.p{wlX + (Xa){xp)} 

= ^dkRk,apWl^X + 9kASk[xl (78) 

where we used the fact that Wk [J] is the generating func- 
tional of connected Green functions. With the aid of H77() 
we immediately obtain a flow equation for the effective 
average action 

dkTklx] = ^dkRk,^p [Ff + Rk] -'m^o. ^^^^ 
= iSTr{afei?fc[F(2^+i?fc]-i}. 

Here the "supertrace" runs over field type, momentum, 
internal indices etc. and has an additional minus sign for 
fermionic entries. 

This equation is exact - we have only performed well 
defined formal manipulations*. However, it is an equa- 
tion for an infinite number of couplings and hence by no 
means accessible to an exact solution. The usefulness of 
eq. H79|) will only show up if we are able to make sensi- 
ble approximations to the flow equation. We will come 
back to this later. Here we note that in the flow equation 
(|79|l the regulator function Rk appears in the effective 
propagator as an infrared (or Fermi surface) regulator 
whereas the factor dkRk reflects the fact that the only 
/c-dependence arises through the regulator. For an ap- 
propriate choice of Rk with dkRk decreasing rapidly for 
large momenta the momentum integrals are ultraviolet 
finite. This means that effectively only a small interval of 
momenta contributes to the integrals. 

Let us rewrite the flow equation in a very useful way 
making contact to perturbation theory. Define the deriva- 



tive (the index i counts the field types) 

~dk = (9kRl)^^. (80) 

With the aid of this derivative the flow equation can be 
cast into the form 

dkTkix] = iSTr{4ln[Fi'^ +Rk]}. (81) 

This can to be compared with the regularised perturba- 
tive one loop result ^ 

TfcM =rA[x] + iSTr(ln[5(2)+i?,]-ln[5(2)+i?A]). (82) 

Performing the fc-derivative of eq. H82(l leads to a one 
loop flow equation. A "renormalisation group improve- 
ment" iS'*-^-' F^^'' promotes this equation to a non- 
perturbative exact flow equation ^20]. This comparison 
with perturbation theory allows us to identify the right 
hand side of H81|) as a sum of one particle irreducible one 
loop diagrams, where all couplings have been replaced 
by their renormalised counterparts. Momentum integra- 
tions, sums over internal indices etc. are performed after 
the dk derivative. 

Obtaining the flow equation for some coupling thus 
amounts to summing all one loop diagrams for this cou- 
pling, evaluating the dk derivative and then calculating 
the trace. However, we may be able to perform parts of 
the trace first if the cutoff does not depend on it. For ex- 
ample we will later be able to first sum over Matsubara 
indices before performing the dk derivative. Technically, 
this is a very useful property since it allows us to use 
directly the one loop computations (e. g. sect. VIII) for a 
derivation of the appropriate renormalization flow equa- 
tions. 

On the level of the exact flow equation we have, of 
course, to remember that the one loop diagrams have 
to be evaluated for field dependent propagators as given 
by (r^^^ + Rk)~^- Therefore the flow equation lO is 
a complex differential equation for functionals. In this 
paper we tackle it by expanding the effective action in 
powers of the fields 

r'=M = E;^Exo.---Xa„ri;l^,...,„„. (83) 

n— ai 

The flow equations of the n-point functions F^"^ can eas- 
ily be derived from eq. H79|l by appropriate functional 
derivatives. However, the flow of some n-point function 
will in general contain higher n-point functions. This is 
a general feature: if we perform a systematic expansion 
of the effective average action, the set of flow equations 
will not be closed. We have to truncate the expansion at 
some point. 



If needed, one may work in a finite volume such that the mo- 
mentum sum is discrete. 



Since we have regularised the propagators we may send A to very 
large values. 
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X. REBOSONISATION OF FOUR FERMION 
INTERACTIONS 

Any partially bosonised theory will generate four 
fermion interaction terms under a renormalisation group 
step corresponding to the diagrams in Fig. 5. 



\ / 
\ / 

A 
/ \ 
-L i_ 



FIG. 5; Box diagrams for four fermion interactions. Solid lines 
denote fermions, dashed lines bosons. 

However, we wanted to capture the complicated be- 
haviour of higher fermion vertices in the bosonic lan- 
guage - this is what the bosonisation procedure is all 
about. One might suspect that it should be possible to 
rebosonise the fermionic coupling obtained after some 
renormalisation group step by a suitable field redefini- 
tion of the bosonic fields. This is indeed the case as was 
shown in (see also fT^I. 

Consider a theory with (effective average) action 



Q 

+ i^(/.(-Q)F^,fc(g)</)(Q) 
Q 

-J2hk{Q)HQ)k~Q) 

Q 

+ ^A^,fc(g)<^(Q)0(-g), 



(84) 



where is the fermionic bilinear corresponding to the 
bosonic field 0, e.g. 0j(Q) = ^j^^p* {K)(t''^:{Q + K), and 
the initial condition for the purely fermionic coupling is 
= at some initial scale k. 
Now perform a renormalisation group step from the 
scale k to the scale k — k — Ak. The change in scale, Afc, 
is supposed to be so small that the changes in couplings 
are also small; they are calculated by the flow equation 
()79|l for the truncation H84(l . As can be seen from the 
diagrams in Fig. 5, the four fermion coupling A^^/c will 
in general be different from zero, say AA^^fc, at the new 
scale k. 

We will use our freedom in the definition of the bosonic 
fields to consider a field redefinition at the scale k 



MQ) ^ MQ) + ^MQ)m), 



(85) 



where Aak is an up to now arbitrary function. (We put 
(j>i: = (f> at the initial scale.) For infinitesimally small shifts 
this results in a flow equation for a ^-dependent fleld 
variable (pk 



Let us see how we can implement the field redefinition 
into the renormalisation group formalism. In the exact 
renormalisation group equation H79|l the change of scale 
9ferfc[x]|^ is calculated at fixed fields. Hence, if in addi- 
tion we perform a shift in the fields as above and consider 
the flow for flxed (pk, the flow equation becomes 

+ E(3lrr[V',^*,</'fe])5fe0fe 

Q 

=dkTk[ij,r,MW (87) 

+ dkakiQ)P4,,k{Q)MQ)k~Q) 
Q 

+ hk{Q)dkMQ)HQ)k~Q))- 

It is appropriate to define the couplings etc. as the 
coefficients of polynomials in 0^. The secons term in eq. 
()87|l this changes the flow equations for hk and X^^k to 



dkhkiQ) = dkhkiQ)\4,^ + dkak{Q)P4,^k{Q), 
dkX^,k{Q) = dk\,j,^k{Q)Uu + hk{Q) dkUkiQ). 



We may now demand that the purely fermionic coupling 
vanishes for all scales fc, i.e. dkX^,k{Q) = 0. This deter- 
mines dkOikiQ) and leads to a modified flow equation for 
the Yukawa coupling 



dkhk{Q) - dkhk{Q)Uu - ^T^^^kX^AQ) 

n-kW) 



(89) 



As a net result, we have traded the "regeneration" of a 
four fermion coupling against a modification of a running 
Yukawa coupling. 

Even though no direct four fermion coupling is now 
present in the truncation H84|l (A.0^fc = 0), the effective 
four fermion coupling can always be reconstructed from 
the partially bosonised effective action. For this purpose 
one solves the field equation for ipk as a functional of ip 
and reinserts the solution ip^^^ [ip] into the effective action. 
For the truncation (|84ll this yields 



4"' 

\eff 



{Q)^hk{Q)P^',{Q)m, 



(90) 



dkCpk{Q)^~dkak[Q)i[Q). 



(86) 



Adding the various bosonic channels which have been 
taken into account in the bosonisation procedure a suit- 
able renormalised local quartic fermion coupling Uk can 
be defined as a linear combination of X'^kiQ ~ 0). Once 
the flow of hk and P^^^ is computed one can compare the 
flow of Uk with the flow coniputed directly in a purely 
fermionic formulation ^ We have checked that 
with our rebosonisation prescription the "fermionic" and 
"bosonic" flow of Uk agree in lowest nontrivial order (e.g. 
dkUk = cC^D- We emphasize that this "one loop accu- 
racy" of the flow of Uk would not hold had we omitted 
the "rebosonisation correction" H89|l . 
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Beyond one loop order the "fermionic" and "bosonic" 
flow of Uk shows differences that depend on the pre- 
cise truncation. (Without truncations the two versions of 
the flow have, of course, to coincide precisely.) They re- 
flect the emphasis on different aspects of the model. The 
bosonic flow investigated here takes also into account the 
effect of the bosonic fluctuations. They contribute to the 
generation of quartic bosonic couplings which would cor- 
respond to eight-fermion-vertices in a purely fermionic 
language. These quartic bosonic interactions are crucial 
for our ability to explore spontaneous symmetry breaking 
by the renormalisation group flow. Indeed, they guaran- 
tee that the bosonic effective potential remains bounded 
from below even if the bosonic propagator P(f,.k{Q) turns 
negative for some range of Q. From equation H9()|l we learn 
that the effective four fermion coupling Uk will diverge 
precisely when P^.fc vanishes. This divergence limits the 
the effective four fermion coupling of the purely fermionic 
flow to the symmetric phase. Going beyond, and in par- 
ticular exploring the ordered phase, would require to take 
at least eight- fermion vertices into account, rendering a 
purely fermionic description complicated. These difhcul- 
ties are overcome by the bosonic flow equations which 
can easily follow the flow of a four-boson-vertex. 

On the other hand, we observe that the four fermion 
vertex multiplied by A^fj. (eq. H9U|) ') shows a particu- 
lar structure in momentum space: from the combination 
X]q ^^,k{Q)4'{Q)4'{^Q) the momentum of a pair of two 
fermions sum up to a common "pair momentum" Q. This 
"factorised structure" does not exhibit the most gen- 
eral momentum dependence of the effective four fermion 
vertex permitted by the symmetries. It is therefore not 
possible to rebosonise the most general momentum de- 
pendent four fermion interaction generated by the flow. 
As long as one uses only local bosonic fields only the 
"dominant part" of the four fermion interaction - pre- 
cisely the one that can be described by boson exchange 
in several channels - can be made to vanish for all k 
with the help of field redefinitions. In our approxima- 
tion, we simply omit the "subdominant part" that does 
not take a factorised form. This contrasts with the work 
on the fermionic renormalisation flow [H Qi |M |^ where, 
in principle, the full momentum dependence of the four 
fermion interaction is taken into account. (Of course, one 
could combine the advantages of both approaches by ex- 
tending the truncation, keeping explicitly the influence of 
the "subdominant" four fermion interactions.) In short, 
whereas the fermionic formulation typically allows for a 
higher resolution of the momentum dependence of the 
four fermion vertex, the partially bosonised version per- 
mits the inclusion of higher vertices and the investigation 
of spontaneous symmetry breaking. 



XI. RENORMALISATION FLOW OF THE 
HUBBARD MODEL 

We now want to apply the renormalisation group for- 
malism summarized in the preceeding sections to the 
Hubbard model. We will concentrate on the region close 
to half filling and low temperatures where the system is 
dominated by the antiferromagnetic spin density. In this 
section we will therefore leave aside all bosons apart from 
the spin density m(X). There is then no ambiguity how 
the parameter Um is related to the original fermionic cou- 
pling U, namely Um = U/3. We extend this truncation 
by including other bosonic degrees of freedom in section 
IX VII We recall this in this setting the mean field approx- 
imation is not consistent with Hartree-Fock or lowest or- 
der perturbation theory {Um 7^ U/2). The consistency 
will be established only by the inclusion of the bosonic 
fiuctuation effects. 

Truncation 

Let us now try to define a suitable truncation for the 
effective action. The initial condition (at very large k) 
for the fiow equation H79|) is the classical action. In the 
course of the fiow towards lower scales the effective aver- 
age action will in general pick up all possible couplings 
that are compatible with the symmetries of the theory. 
In order to make progress we have to truncate this set of 
infinitely many couplings. We will make an ansatz con- 
taining a fermionic kinetic term F^^^, a term containing 
a Yukawa like interaction between fermions and bosons 
Ty^k and a bosonic action to be specified below. (A term 
containing a four fermion interaction is to be rebosonised 
as sketched in the previous section.) As we are mainly 
interested in antiferromagnetic behaviour we define the 
boson 

a(Q)=TO(Q + n), (91) 

whose zero momentum mode a(0) corresponds to a ho- 
mogeneous antiferromagnetic spin density. 

For the fermionic kinetic term we adopt the classical 
part unchanged 

F^.fciV-,^*] =^V^(Q)Pf(Q)V'(Q), (92) 
Q 

Pf{Q) — i^F + e — /X, e(q) = ~2t{cosqx + cos^j^), 

where we restrict ourselves to nearest neighbour hopping. 
In contrast, the Yukawa coupling ha,k is taken to be scale 
dependent 

Fy,fe[V',V*,a] = -Vfe II d{K)r{Q)S^{Q') (93) 

KQQ' 

x5{K -Q + Q' + 

For the purely bosonic part we investigate a kinetic term 
plus a local effective potential 

Q X 
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Due to SU{2) symmetry the potential can only depend on 
the rotation invariant combination a{X) — ^a{X)a{X) 
or, in momentum space 



a{K,K') = ^a{K)a{K'). 



(95) 



Furthermore we will make an expansion in powers of the 
field a up to a quartic interaction. We take a different 
expansion point in the symmetric regime (SYM) and in 
the regime with spontaneous symmetry breaking (SSB) 
as it is preferable to always expand around the minimum 
of the potential 



SYM : J2 = ^m2a(_i^,i^) + (96) 

X K 

+ 1 '>^aSiKi + K2 + Ks + K4) 

Ki...Ki 

SSB:^L/[a] = \ ^aSiKi + K2 + K3 + K^) 

X K1...K4. 

x{a{Ki,K2)-aoS{Ki)S{K2)) 

The parameters in the potential are the scale dependent 
mass rria{k) or the minimum ao(fc), as well as the cou- 
pling Aa(fc). For large k we start in the SYM-regime and 
switch to the flow equations for the SSB-regime at the 
scale fcsSB defined by m^(A:ssB) = 0. 

The bosonic propagator on the classical level is simply 
a mass term in our case, i.e. Pa — Xa = 0. We let us 
guide by the one loop results of the preceding section for 
the momentum dependence and take 



(97) 



where Za is a scale dependent wave function renormalisa- 
tion (with the dimension of mass^^) and the function [qY 
is defined as [q]^ — ll + l^ for qi G [— tt, tt] and continued 
periodically otherwise. 

The ansatz [q]^ for the spatial part is to mimic the 
momentum dependence of figure|31 It is not very accurate 
for 5? in the vicinity of tt^. For k A this region gives 
only a small contribution, however, whereas for k near 
A the whole kinetic term is small as compared to the 
mass term. (The classical approximation is Za = 0.) The 
dependence on the Matsubara frequency is more difficult. 
A look at figure 0] shows that the bosonic propagator 
is changed only for the smallest frequencies, while the 
higher frequency modes retain the classical mass term 
(without sizeable kinetic term) . This could be reflected by 
the ws-dependent part Xa{ujB) m Pa obeying Xa{0) = 0, 
Xa « m^j — m^(fc) for large ujb- Lowering the scale, the 
ujb ~ mode will then dominate the propagation more 
and more. Our approximation Xa{iOB) = ZaijJ% chosen 
for the sake of simplicity: it has the correct qualitative 
property of suppressing the effect of the high frequency 
modes. 



One would suspect a similar behaviour for the quartic 
boson coupling Aq which is generated by the fermion loop: 
the low frequency modes are supposedly changed most, 
implying that Xa depends substantially on the frequency 
lob- (For example we could replace A^ 11^=1 /a('j-'b) 
with /a(0) = 1 while f\ — for large usg, where i labels 
the external legs of the diagram contributing to the flow 
of A.) We neglect this effect here since the contribution of 
modes with large wb is suppressed anyhow by the term 
^ ZaLu% in the inverse propagator. 

In addition to the truncation we still have to spec- 
ify the regulator functions for the renormalisation group 
equations. 

Fermionic regulator 

The fermionic cutoff function is inspired by the 
fact that at nonvanishing temperature the propagator 
Pf{Q) — i^F + e — n has no zero-modes. This means 
that the temperature itself acts as a regulator. We there- 
fore choose 

RtiQ) = iuFi^ - 1) = 27:iinF + i)(rfe - T), (98) 

which has the effect of replacing the temperature T by 
some function Tk in the fermionic propagator. We specify 
this function to be {p > 2) 



rpP _ rpp 



kP 



dkTk 



p-i 



\{k/T)P-^ iffc 



if fc>r 



(99) 
(100) 



such that for large k the flow mimics a change in temper- 
ature whereas for k <^ T the modification of the fermion 
propagator becomes ineffective. For fermionic modes for 
which the distance of the momentum to the Fermi surface 
is smaller than ttT no additional regularisation is needed 
anyhow. 

In short, we want a cutoff that is ineffective for fc <C T 
and allows us to approach the Fermi surface continuously 
for k > T. We will see in the next section that this re- 
quires p > 2 and we will choose p = 4. 

Bosonic regulator 

For k <^T possible infrared problems may arise from 
the fluctuations of long range bosonic modes. This is par- 
ticularly relevant near a second order phase transition 
where the bosonic correlation length diverges. Such po- 
tential problems are avoided by the bosonic regulator. 
We take a "linear cutoff" [2I 



i?^(g) = Z,(fc2-Q2)0(fc2 



(101) 



where Cp is defined in eq. H97|l . This leads to a full prop- 
agator of the form 

PkiQ) = Pa{Q) + Rl{Q) (102) 

where 



Ql - 



Q2 if Q2 > ^2 



(103) 
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The regulator function (|101|l thus only hampers the prop- 
agation of modes with small momenta and Matsubara 
frequencies where < k^. Therefore, by lowering k, we 
effectively average over larger and larger regions in po- 
sition space. We may therefore relate properties of the 
effective average action Tk at a given scale k to prop- 
erties of size 1/k in position space. Since the fermionic 
modes are effectively integrated out for k <^ T we can as- 
sociate in this range k~^ with a typical (linear) size of the 
experimental probe. Indeed, fluctuations with size larger 
than the experimental probe should not be included in a 
"physical" definition of the relevant effective action. For 
practical purposes one is therefore interested in macro- 
scopically large but not infinite values of A;^^. 

Initial values 

For k oo both the fermionic and bosonic propaga- 
tors vanish. Therefore the loop contributions vanish and 
rfc^oo indeed equals the classical action. The initial con- 
ditions for the parameters occurring in the truncation 
can therefore be read off from the classical action (|43|1 . 
For very large k one starts in the SYM regime with 



U„ 



Xn 



0, ha 



Ur, 







(104) 



XII. FLOW EQUATIONS AT HALF FILLING 

In this and the following section we study the Hubbard 
model at half filling, i.e. for vanishing chemical potential 
/i = 0. As a test of the reliability of our approach we ex- 
tend the truncation in sect. IXVU and check the influence 
of the Fierz ambiguity. 

Bosonic potential 

The flow equation for the effective potential U{a) ob- 
tains by evaluating the ERGE (|79|l for a homogeneous an- 
tiferromagnetic order parameter d{Q) = ad{Q), a = ^a^. 
The contribution of the fermionic fluctuations can be 
found by replacing in the mean field theory result H53() 
T — > Tfe in the fermionic propagator and applying the 

-derivative. The bosonic contribution is the same as 
for the 0(3) Hnear cr-model 



dkU{a)^dkU''(a) + dkU^{a) 

= i 5] 9fe ln[F,(g) + Mfia) + RHQ)] 



(105) 



2T 



:dk In cosh y (a). 



Here the boson mass terms are different for the "ra- 
dial mode" (in the direction of a) and the "Goldstone 
modes" (perpendicular to a), 

+ 3Aaa, ml + Xa, m^-f !„«) SYM 
\ (Aa(3Q! - ao), Aa(a - ao), Xaia - ao)) SSB 

(106) 



In the SSB regime we recognize the two massless Gold- 
stone bosons for a = Uq. The squared fermion mass (mass 



gap) 2/i„a enters the function 



1 

m 



2hla. 



(107) 



The flow of the mass term and quartic coupling follow 
from simple differentiation of eq. (|l()5|l 



SYM: dkmi 



i:{dkU{a))U=o, 



SSB: 



dkXa = ^(<9feC/(a))|a=o, 
dkaa = -^-i:{dkU{a))\c=ao, 



(108) 



dkXa = ^{dkU{a))\ 



Here the flow of the minimum is inferred from the fact 
that U'{ao) = and hence ^t/'(ao) = dkU'{ao) + 
U"{ao)dkao = 0. 

With the aid of 1)1081) we may now derive the flow 
equations for the parameters in the effective potential. 
However, first we introduce rescaled and renormalised 
quantities 



such that 



Zak^ 



Tk^' 



du 
da ' 



T 



(109) 



T 



T - _ d^u 
zWk^ ~ d^' 



f,2 T -2 



(110) 



Zat 



4 a 



are all dimensionless. 




FIG. 6: Flow of the renormalised and unrenormalised mass 
10~^m^ (solid), lQ^m1/t (long dashes) and the quartic 
bosonic coupling \a (short dashes) in the symmetric regime 
(SYM) at half fiUing. We have chosen U /t = 3 and T /t = 0.15. 



The flow equation for u can be expressed in terms of 

_2 

"threshold functions" {w) which express the decoupling 
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of particles with mass much larger than k (large w) 

kdku\a = -2u + r]aau' + j^kdkU\a 

= - 2u + rjaau' + ^^ll{wi,r]a,T) 



+ So{rhp,T) 



(111) 



This uses dku\a — dku\a + rjadtu' , u' — du/da and ija is 
the anomalous dimension 

r]a = -kdk\nZa. (112) 

The dimensionless threshold functions (g — tq) 



s^{^%,T)=[-] (- 



(113) 



1 /2 

y = (4(cos9i + cos 92)^ + m%) 

decay rapidly for large values of the arguments w — 
lVP/{Zak'^), fri\ = 2hla. Their explicit form IS given m 
the appendix. The characteristic scale dependence of the 
mass term and quartic coupling is shown in Fig. El The 
running for k ^ t has not much physical meaning and 
rather sets the stage for the correct translation between 
the classical action and F^-f . The physical fluctuation ef- 
fects become important for k < t as seen in the running 
of m^. ^ ^ 

It is instructive to consider the long wavelength limit 
fc^ -C T^. For p > 2 the fermionic contribution Sq van- 
ishes ^ kP^'^. As argued in the previous section the 
fermionic infrared cutoff should be ineffective for k <^T . 
This means that the change in becomes irrelevant. 
In turn, this requires sufficiently large p and we choose 
p = 4. On the other hand the bosonic threshold function 
approaches for k <^ T the one for the two dimensional 
linear 0(3) cr-model HHlEi 



lQ{w,ria,T) lliw,ria) 



1 - Va/i 

1 + w 



(114) 



In the long wavelength limit we therefore recover the well 
studied flow equations for the scalar 0(3) model in two 
dimensions |2q| . The flow equations for rri^, and k 
involve the a-derivatives of Iq and Sq and can be found 
in the appendix. We note that all these equations have 
a simple diagrammatic representation. For example, the 
mass contributions in the symmetric regime are shown in 
Fig.CI 

Anomalous dimension 

The anomalous dimension 77^ = —kd^la Za is a mea- 
sure for the change of the wave function renormalisation 




FIG. 7: Diagrams for the running of the scalar two point func- 
tion. 



Za with scale. Therefore we can extract it from the mo- 
mentum dependence of the bosonic two point function. 



Ts.kin = Y,Pa{K)a{-K,K) = Z K'' a{- K , K) . 

(115) 



K 



K 



As we are mainly interested in the spatial momentum 
dependence we define 



d 



t2 d{P) 



(116) 



K = [ljb = 0,k = lei), 



where ao is the minimum of the effective potential. We 
observe that Za is closely related to the momentum de- 
pendence of the inverse propagator shown in figs. 3, 4. 
The definition (|115|l . 111()|) corresponds to an approxi- 
mation near the minimum, where the constant part (at 
the minimum) is substracted. 

In the symmetric phase (SYM) the bosonic propaga- 
tor is affected by the two diagrams in Fig. 7. However, 
the bosonic loop is independent of the external momenta 
and therefore does not contribute to the anomalous di- 
mension. The fermionic loop is well known from our loop 
calculation by the replacements T — > T^, [/„i — > ha for 
Pa{K) = Pm{K + n) in ea. l|M|) . and performing the dk- 
derivative. We obtain 



(F) 

ria 



hltikdkTk) 



a 



tanh ^rr^ -|-tanh 



OTk d(P) 

e(q+iei) 



(117) 



2T,, 



(e(q)+e(q+;ei))/t 



(=0 



A characteristic evolution of Za and rja is plotted in Fig. 

El 

In the S'5i3-regime we also get a contribution from the 
bosonic sector. The contribution comes from a bosonic 
loop with four external legs, where two external legs are 
connected to the condensate. As we will see later, the sys- 
tem enters the broken phase at small values of the cutoff 
parameter k. At these values only the lowest Matsubara 
frequency {ujb — 0) contributes in the bosonic propaga- 
tor (c.f. (|102|l ). The bosonic contribution to rja is then the 
same as the corresponding one for a simple 0(3) model in 
two dimensions. Here the anomalous dimension has been 
calculated to be [201 



n{l + 2\aHaY 



(118) 
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FIG. 8: Flow of the wave function renormalisation Za (solid) 
and the anomalous dimension \(i~^r]a (dashes) in the sym- 
metric regime (SYM) at half filling. We have chosen U /t = ?> 
and T/t = 0.15. 



where the function ra\ 2 {wi , W2 ) contains the momentum 
dependence of the loop integral and depends on the regu- 
lator i?^. For our choice of regulator it has been computed 
exphcitely 01 as shown in ifTll^ . In the SSB -regime the 
contribution from the fermionic loop is much smaller than 
the bosonic contribution and we will neglect it. 

Running Yukawa coupling 

In the symmetric phase the running of the Yukawa 
coupling is generated by the diagrams in Fig. 



the beta functions read {ujp = (2n -I- l)7rTfc) 



id) 

hi 



2ht kdk J2 {Pa ' {Q)pf' {Q)pf\Q + n) } 



■E 



kdkUJp 



1 



^ \ {QjI + e2(q))2 Q2/fc2 + ^2 



1 



(2-r7,(l-Q7fc2))e(fc2_Q2) 



c2;|, + e2(q) 



(l + m2)2 



/Sg"' =2^lhl—m J2 {PaHQ)Pa\^ Q)Pf\-Q) 



1 (2-r;,(l-Q2/fc2))e(fc2-Q2) 



2 1 



kdkUJ-p 

'fc2T^ 1 (t2;2,+e2(g))2Q2/p^^2 



(l + m2)2 

1 



(Q + n)2/fc2 + m2 



(120) 

In cq. inni) we use Ql = Q^Q{Q^ - + PQ{k^ - Q^), 
kdkOJp = (2n + l)iTkPTl^P and 

Pf{Q) = [iLOF + e(q))A, PaiQ) = Ql/k^ + ml (121) 

We plot the flow of the renormalised and unrcnormalised 
Yukawa coupling in the symmetric regime in fig. 10. 
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FIG. 9: Diagrams for the running Yukawa coupling 




-Ink/t 



The first diagram is the direct contribution, while the 
last two contribute by rebosonisation as prescribed in 
sect. 1X1 More details on the extraction of the rebosonised 
flow equations can be found in appendix 

In the symmetric phase we get for the running of ha 

kduhl = hl^, + p\^+l3^l^ (119) 
where the "direct" and "rebosonised" contributions to 



FIG. 10; Flow of the renormalised lQ~^ha (solid) and un- 
rcnormalised ha (dashes) Yukawa couplings in the symmetric 
regime (SYM) at half filling. We have chosen U/t = Z and 
T/t = 0.15. 

For the specific parameters used for this plot we ob- 
serve very little running of ha- The charge will be more 
pronounced for other choices of the parameters, as dis- 
cussed in sect. XV. In the regime with spontaneously bro- 
ken S'J7(2)-symmetry (SSB) the change in the Yukawa 
coupling is negligible as we have checked numerically. We 
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FIG. 11: Flow of the renormalised minimum of the potential 
Ka (sohd), the quartic bosonic coupling 10~^Aa (long dashes) 
and the wave function renormalisation 1Q~^ Za (short dashes) 
in the broken (SSB) regime at half filling. We have chosen 
U/t = 3 and T/t = 0.15. 




10 20 30 40 50 60 70 80 90 

-Infc/t 

FIG. 12: Flow of the renormalised minimum of the potential 
Ka in the broken (SSB) regime at half filling for different tem- 
peratures. We have chosen U/t — 3. The temperatures are in 
the set T/t = {0.05,0.07,0.085,0.10,0.13,0.17} (from top to 
bottom). For T/t = 0.05 the order persists up to fc 10"*"t 
corresponding to a sample size of 10'^^ cm. This clearly demon- 
strates that the Mermin- Wagner theorem is not practically 
applicable. 



neglect the running of the Yukawa-coupling ha{k) in the 
SSB regime. For k < /cssb we therefore keep the unrenor- 
malised Yukawa couphng fixed ha{k < fesss) = ha{kssB)- 



XIII. SPONTANEOUS SYMMETRY BREAKING 

We now turn to a numerical analysis of the above flow 
equations. For this we set U/t = 3 and take a tempera- 
ture T = 0.15^. The initial scale fco — 100^ is chosen so 
large that the final results do not depend on it and the 
one loop results are well produced in the beginning of 



the fiow. The differential equations were integrated by a 
standard Runge-Kutta like routine . In fig. 6 we plot 
the flow of the renormalised mass m^, the unrenormalised 
mass TO^ and the quartic bosonic coupling Aq . The values 
of k correspond to a flow in the symmetric regime (SYM) 
where the minimum of the effective potential occurs for a 
vanishing order parameter a = 0. One observes that the 
running of TOq sets in only for k < t. Then both Wia and 
TUa decrease. At In/c/tw —2.2 the bosonic mass van- 
ishes and for smaller values of k we enter the SSB regime. 
We denote the corresponding scale by fcsSB- In flgure 8 
we plot the flow of the wave function renormalisation 
Za and the corresponding anomalous dimension rja- This 
explains qualitatively the flow of the renormalised cou- 
plings in fig. 6: for scales above Ink/t « 2.0 the running 
is mainly dominated by the simple scaling due to the re- 
spective dimensions of the couplings. In an intermediate 
range down to Ink/t ^ —1 the large value of rja domi- 
nates the flow while for even smaller values of In k /t the 
fermionic part of the flow equations play an ever increas- 
ing role. 

In figure 11 we explore the SSB regime with a{k) ^ 
for k < fcssB- First, we observe that the quartic bosonic 
coupling reaches a fixed point very soon. This is because 
the positive term in equation ljB4|l just compen- 

sates the negative contributions from the fermions and 
dimension-scaling. Two opposite effects govern the flow 
of the location of the minimum of the potential, as param- 
eterized by Ka- For k only somewhat smaller than fcssB 
the fermionic fluctuations dominate the flow. They lead 
to increasing values of Kq. However, soon the fermionic 
contribution becomes smaller due to the increasing mass 
(gap) of the fermions. Then the bosonic loop involving 
the massless Goldstone bosons dominates. This only re- 
sults in a slow logarithmic "running" of Ka , but it finally 
drives the minimum to zero and thus restores the sym- 
metry. When the fermionic part becomes negligible we 
effectively deal with a bosonic 0(3) model in two dimen- 
sions for which the symmetry restoration is a well known 
feature and has been extensively studied in our functional 
renormalisation group formalism 20, 22, 26]. 

This picture reconciles the antiferromagnetic order 
with the Mermin- Wagner theorem , which states that 
a continuous symmetry cannot be broken at nonvanish- 
ing temperature in two dimensions and below. We indeed 
loose the notion of symmetry breaking if we average over 
arbitrarily large volumes, i.e. lower the cutoff parameter 
k to zero. For fc — > the order parameter a{k) always 
vanishes, in accordance with the Mermin- Wagner theo- 
rem. This is due to the presence of Goldstone modes in 
the SSB regime^°. Nevertheless, we find that for low tem- 
peratures the symmetry is restored only when averaging 
over extremely large samples unaccessible to any real ex- 
periment. 



For 0(2) models in two dimensions there is another possibility to 
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For example, for T = O.OSt the restoration of the 0(3)- 
symmetry would require the effects of the Goldstone- 
boson fluctuations with characteristic length scale (in- 
verse momentum) of lO^^cm. For any realistic experi- 
mental length scale, however, the order parameter a{k) 
remains nonzero. This includes macroscopic length scales. 
We conclude that it makes perfect sense to speak about 
long range antiferromagnetic order. 

Let us consider a probe with some finite macroscopic 
size, say 1cm. Obviously, the fiuctuations with momenta 
smaller than (lcm)~^ are absent in this setting and we 
should stop the renormalization flow at some nonzero 
k = {lcm)~^ ~ 10~^t. For all observations with charac- 
teristic length scale smaller than Icm (in order to avoid 
boundary effects from the particular geometry of the 
probe) one observes all features of an effective antifer- 
romagnetic state: For all temperatures smaller than a 
"critical temperature" Tc there will be an average spin 
direction oq ("spontaneous symmetry breaking"). Fur- 
thermore, the fluctuations orthogonal to oq will behave 
as two Goldstone modes with effectively "infinite" (^ 
size of the sample) correlation length for all T < Tc- 
On the other hand, the antiferromagnetic spin wave fluc- 
tuations in the direction of oq (radial mode) will have 
a finite correlation length ^ for T < Tc. As Tc is ap- 
proached from below also the radial correlation length 
diverges ^(T Tc) cx)", i.e. exceeds the size of the 
probe. The different behavior of the Goldstone and ra- 
dial modes is a direct consequence of the presence of a 
direction specified by ap. Another consequence of effec- 
tive antiferromagnetic order is the presence of a mass gap 
for the electrons and holes. 

Let us next turn to the notion of a critical tempera- 
ture Tc. We will see that its precise value depends on the 
size A;"^ of the experimental probe. For this purpose we 
display in Fig. 1121 the running of Ka for various temper- 
atures. For a fixed value of fc > (corresponding to the 
inverse of the sample size) we always observe a transition 
to an ordered phase if T is low enough. This transition 
temperature is denoted by Tc{k). In fig. 1 it is directly 
visible as the temperature for which the order parameter 
Ka vanishes. In fact, fig.Qshows Ka{k) with k/t — 10~^ 
much smaller than all other characteristic scales of the 
model. 

Inspecting fig. E| we may define a characteristic scale 
oof symmetry restoration fcsR by the vanishing of Ka, 
i.e. (a(fcsR) = 0. Only "domains" with a size of kg^ be- 
come so randomly oriented that no net order prevails on 
length scales larger than kg^ . For an experimental setting 
with characteristic scale k ^ ksn the formal "asymptotic 
symmetry restoration" for fc is without practical rel- 
evance. Thus, even though formally correct, the Mermin- 



Wagner theorem fails to be practically applicable in the 
case of the Hubbard model at low temperature. In con- 
trast, for fc ^ fcsR the precise value of fc is unimportant 
and we may take fc = 0. In fact, this situation corre- 
sponds to the ordered phase with a renormalized mass (or 
inverse correlation length) mR ^ fcsR ^2,6]. The correla- 
tion length ^ acts as an effective physical cutoff such that 
fluctuations on length scales larger than ^ are irrelelvant 
and an additional "experimental cutoff" k <^ mji has no 
effect. 

XIV. CRITICAL BEHAVIOR 

For T sufficiently large compared to Tc our analysis 
confirms the statement of that the long distance be- 
havior of the Hubbard model can be described by a clas- 
sical nonlinear cr-model. This holds for T in the vicinity 
of and below Tc. We can use our findings in order to 
establish a quantitative description of the temperature 
dependence of the correlation length ^ ~ mjf^ . For char- 
acteristic momenta ^ T the fermion fluctuations are cut 
off by the temperature and only the bosonic flutuations 
remain relevant. Within the remaining effective classical 
0(3) linear cr-model the running of k obeys for k suffi- 
ciently large [26l | 

fcSfeK- -^ + — V+0(«-2). (122) 

As a reasonable approximation (cf. fig. Ill|l we may em- 
ploy a "linear" behavior for fc between fc™ and fcsR 

K(fc) = At™(r)--^ln^^ (123) 
47r fc 

where (T) corresponds to the maximum of the curves 
in fig. ^ With = fcsR/CsR and K(fcsR) = one 
obtains for the correlation length 

e = £^exp(4^Ac„(r)). (124) 

Using the definition UllOfl for k this yields the character- 
istic behavior js^l 

e = Cexp(^) (125) 

with 

J ^ AnaoiT)Za{T)t\ (126) 

Here cS-^^Z^^ denotes the T-dependent maximum of 
(ao^a)(fc) and both 7 and C = Csn/km depend smoothly 
on T (without particular features for T — s- 0). 



circumvent the Mermin-Wagner theorem related to a KosterUtz 

and Thouless phase transition fSl. It is speculated that this kind " The nonlinearities result in a proportionality factor fesR = 

of mechanism may play a role in the superconducting region H fesRCsR and in a slight derivation of km from the location of 

1221 • the maximum in fig. 1121 
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For T < Tc the correlation length according to eq. 
(|125ll exceeds the size of the experimental probe and eq. 
(|125|l is no longer applicable. The modifications of the 
results of for T < Tc therefore concern the tem- 
perature region where the formal correlation length (for 
fc ^ 0) would become extremely large and not practically 
meaningful anymore. We can actually use eq. H123|) for 
an estimate of the critical temperature Tc{k) 



Tc{k) = 



liTc) 



In (k,n{Tc)/k) 



(127) 



The fc-dependence of Tc is only mild - typically a change 
in fc by a factor of 10 changes Tc by 10 percent. For a sam- 
ple size k = 10~^i ~ (1cto)~^ we may use a quantitative 
formula for ^ 



= c(T) exp <^ 20.7/3(r) 



with 



c{T) 
5{T) 



T 

c 

Y 



ao{T)Za{T) 
ao{Tc)Za{Tc) ' 



(128) 
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km{T) 

Tr 



5{T) 



T 



1 



(129) 



and P{Tc) = l,c(Tc) = Csr. Despite the fact that Tc 
is not a sharp temperature the behavior for T near Tc 
shares several features of standard critical behavior. It is 
obvious that the correlation length increases very rapidly 
as T is lowered from say 1.5Tc to l.lTc - still, it does not 
diverge for T Tc- On the other hand, for T < Tc 
the order parameter reaches zero for T Tc (fixed k) 
according to 



«a(T) 



liT)Tc 
T 



1 



ATc 



1 , kr. 

— in — 

47r k, 



.{Tc) 



(130) 



The critical temperature Tc — Q.llAt can be directly read 
off from fig.n Furthermore, in the vicinity of Tc the prop- 
agator G{q^) for the antiferromagnetic spin waves with 
momenta 3> i^^^ is characterized by an anomalous 
dimension rja — rja{k = \fcp-^T = Tc), namely 



flow enters the SSB-regime. We will call this the pseu- 
docritical temperature Tp^ and infer T^c = 0.18i (for 
V jt — 3). For T > Tpc the mass term m1{k) remains 
positive for all k. We note that Tpc is substantially above 
Tc {Tpc/Tc ~ 1.6). In the temperature range between Tc 
and Tpc we may describe the physics by antiferromagnetic 
domains with fluctuating directions. The typical size of 
domains with local order extends up to the scale kg^ 
where Ka{ks^) = 0. On length scales larger than fcgj^ no 
long range order remains. It is obvious that the critical 
behaviour described by eq. (|123|) looses its validity for T 
above Tpc or in the vicinity of Tpc- Therefore Tpc roughly 
describes the upper bound of the range of validity in T 
for the critical behaviour discussed above. 

The positivity of for T > Tpc implies that the effec- 
tive four-fermi-interaction mediated by the exchange of 
the a-boson remains flnite. We conclude that the "criti- 
cal temperature" determined by the fermionic flow equa- 
tions 3, JL, _5, J)] actually corresponds to Tpc rather than 
Tc- Indeed, the approach based on the fermionic flow 
equation precisely looks for the "critical temperature" 
where the four-fcrmion-interaction diverges for some k - 
or, in our language, where reaches zero. The same 
holds for mean fleld theory or the Hartree-Fock approx- 
imation. Again, these approaches determine Tpc rather 
than Tc- The true critical behaviour discussed in this 
section is governed by the renormalisation flow for the 
effective interactions of composite Goldstone bosons. It 
cannot be captured by the Hartree-Fock approximations 
or the fermionic flow equations. We emphasize that the 
influence of the Goldstone boson fluctuation on the value 
of the critical temperature is a large effect - we find Tc 
almost a factor of two lower than for the mean field 
theory. 



XV. RENORMALISATION FLOW FOR 
NONZERO DOPING 

In this section we want to compute parts of the phase 
diagram for nonzero doping. For this purpose we study 
the flow equation for nonvanishing chemical potential in 
the symmetric phase (i.e. /i ^ 0, a = 0). The bosonic part 
of the flow of the effective potential is not altered, while 
the fermionic part gets modified. Similar to eq. (|lll|l we 
obtain 



G(g2) = (Z,(fc = ^q^)q 



2^-l+r;„/2 



(131) kdku\a 



We emphasize that all the critical features discussed in 
this section extend beyond the particular approximations 
of this work and beyond the particular next-neighbour- 
coupling Hubbard model. They are universal and hold 
whenever the long distance behaviour can be described 
by an effective two dimensional model for a spin field 
with 0(3) symmetry. 

Inspection of fig.ll2lreveals another characteristic tem- 
perature, namely the minimal temperature for which the 



2u + rjau' 
^5I^o(w^,f?a,T) 



(132) 



A substan tial lowering of Tc has also been observed in a 1/d 
expansion l31l . with d the dimensionality of space. The precise 
connection with our result remains unclear, however, since the 
1/d expansion at most partially captures the Goldstone boson 
physics. 
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where Iq is defined as in eq. while Sq reflects the 

different location of the Fermi surface for /i 

^''?,9jlncosh^ + lncosh^}, 



k (2^)2 - ^ 2Tk 

1 /2 

y = (4(cos(7i + cosga)^ + "T-l) • 



(133) 



The flow equations for the couplings and Aq in the 
bosonic potential can again be obtained by appropriate 
derivatives with respect to d (cf. eq. (|110|l '). 

The equation for the anomalous dimension becomes 



rja = hltkidkTk) 



_d d_ 

(pq tanh 



(134) 



2Tk 



tanh 



e(g+iei)-jj 
2Tfc 



{e{q) + e{q + lei))/t 



1=0 



For the Yukawa-coupling we obtain the same functional 
form of the beta functions as in eq. p2U|) . However, the 
fermionic propagator is now replaced by pf{Q) = {iijJr + 
e(q) -At)/t (cf. eq. 

We have analysed the phase diagram of the Hubbard 
model for small values of the chemical potential ^. The 
results are plotted in figure 13, where again we have cho- 
sen U — 3t. The upper line shows the temperature Tpc 
at which the bosonic mass vanishes in mean field the- 
ory (cf. eq. (|64|l for Q — H). For the bosonisation we 
have chosen Up — 0. For the small values of /i considered 
here the MFT approximation yields a second order phase 
transition such that the critical line indeed corresponds 
to a vanishing bosonic mass term. (In order to deal with 
first order phase transitions one would have to treat the 
bosonic potential in a more complicated truncation. We 
therefore restrict ourselves to small values of // here.) In 
the MFT- approximation Tpc(/i) therefore indicates the 
phase boundary and Tpc coincides with Tc- 

The lower curve in figure 1131 shows the pseudocritical 
temperature Tpc for various values of the chemical poten- 
tial derived with the aid of the flow equations displayed 
above (again with U — it, Up — 0). One observes that 
the pseudocritical temperature is lowered as compared 
to the mean field result. This is precisely the effect of 
the bosonic fluctuations. Indeed, the bosonic fluctuations 
have the tendency to stabilize the symmetric phase - they 
give a positive contribution to the mass term. This can 
explain why the pseudocritical temperature is lowered. 
Another effect is the influence of the bosonic fluctuations 
on the running of the Yukawa coupling. The larger the 
Yukawa coupling the stronger is the destabilisation of the 
symmetric phase by the fermion fluctuations. This effect 
can lower or enhance Tpc, depending on the precise choice 
of the bosonisation. It will be discussed in more detail in 
sect. IXVII The effect of the bosonic corrections is only 
moderate 10%) for the parameters of our example. 
This can be very different for other choices of partial 



FIG. 13: Pseudocritical temperature Tpc versus chemical po- 
tential /X for = 3 in the mean field approximation (above) 
and with flow equations (below). 



bosonisation, as visible in fig.^jfor large Up/t where Tpc 
is lowered by almost a factor of 2. 



XVI. REMOVING THE MEAN FIELD 
AMBIGUITY 

In this chapter we want to investigate how well the 
inclusion of running couplings is able to solve the am- 
biguity with respect to the choice of Yukawa couplings 
in the bosonisation procedure, which was so annoying in 
the mean field calculation. For this purpose we add in our 
truncation the p-boson corresponding to fluctuations in 
the charge density. To keep things simple, however, we 
reduce the bosonic effective potential to a simple mass 
term for each boson and restrict our discussion to the 
symmetric phase. 

For the full inverse propagator P^{Q) = Pa{Q) + 
RliQ) of the a boson we choose 



Za[qf 



for 7^ 
P^{q) fortJB-0 



(135) 



and similarly for P^{Q) (however, we fix Zp = 1). The 
function [q]^ is defined as below eq. (|97|l . This choice re- 
fiects the fact that in the one loop calculation we found 
that the lob = mode is changed most. Furthermore, if 
we make an ws-independent choice of the bosonic reg- 
ulator, we are able to perform the Matsubara sums in 
the loops for the Yukawa-couplings, which drastically 
speeds up the numerics. We choose here a masslike cutoff 
— for both a and p. The fermionic kinetic part of 



the truncation is chosen as in section IXH Specifically, we 
restrict ourselves to nearest neighbour hopping. Further- 
more we will only consider /i = 0. Also the Yukawa part 
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is chosen as in section [XFl Details of the flow equations 
for this truncation can be found in appendix IB 31 

This truncation is a very primitive one. It was chosen 
in order to highlight the dominant factor for the cure of 
the Fierz ambiguity, namely the inclusion of the bosonic 
fluctuations beyond mean field theory. The simple choice 
also makes the numerics relatively fast. We therefore do 
not expect very precise results but one should be able to 
see the general features of the flow. 

The first check was to see if the one loop corrections to 
the fermionic coupling U are reproduced by the above 
flow equations. For small U , and T sufficiently above 
Tc, the perturbative expansion in U is valid and one 
may compute the correction Af7 ~ cU^ in the purely 
fermionic language. Here we compute the same quantity 
in the partially bosonized version. For this purpose we 
start at a large value of the cutoff k, follow the fiow for 
a while and "integrate out" the bosons with their new 
couplings to obtain the new U. As mentioned before, 
this amounts to a solution of the scalar field equations in 
presence of fermionic "background fields" . As expected, 
it turns out that for large enough values of k the flow 
indeed obeys the correct one loop equation. For small U 
and large T we recover the perturbative result. This holds 
independently of the choice of the bosonisation param- 
eters Ui. This is of no surprise because of the one loop 
form of the flow equations and the fact that our trunca- 
tion and rebosonisation systematically includes all terms 
contributing in order U"^. We emphasize, however, that 
the rebosonisation of sect. His needed in order to obtain 
one loop accuracy |19j. 

If we follow the flow towards smaller values of k, the 
purely fermionic coupling will not remain local but rather 
obtain a complicated momentum dependence. Further- 
more, the loop calculations are no longer adequate as a 
comparison for the quality of the flow. We therefore need 
another quantity to investigate the invariance of the flow 
under different parameterisations of the bosonisation. For 
this purpose we have chosen the pseudocritical tempera- 
ture Tpc where the mass of the boson a vanishes. 

In figurelTHwe plot the pseudocritical temperature Tpc 
for the transition to the antiferromagnetic phase for dif- 
ferent values of the parameters Um and Up. Thereby we 
keep a fixed value of the four fermion coupling U = 12t = 
3Um — Up. This graph has been calculated both with the 
flow equations and in the "mean held approximation" , 
i.e. by searching for the zeroes of the bosonic mass in the 
fermionic one loop calculation of the two point function 
The fermionic cutoff was chosen to be — T'^ + fc^ 
for this plot. One observes that even for the solution of 
the flow equation the pseudocritical temperature still de- 
pends on the choice of the bosonisation. However, this 
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FIG. 14: Pseudocritical temperature Tpc for different choices 
of the parameter Up, calculated with flow equations (solid 
line) and in the mean field approximation (dotted line). The 
fermionic coupling is U = 12t = 3Um — Up. The Hartree-Fock 
(or Schwinger-Dyson) result corresponds to the mean field 
value at Up = 6t, i.e. Tpc « 2.9t. 



dependence is relatively mild as compared to the mean 
field result and certainly due to our poor truncation. This 
shows that even a crude inclusion of the bosonic fiuctua- 
tions leads to a substantial improvement over mean field 
theory as far as the Fierz ambiguity (choice of Ui) is con- 
cerned. The comparison in fig. with the Hartree-Fock 
result (corresponding to Up = Um — 6t) also demon- 
strates the overall size of the change in Tpc due to the 
bosonic fluctuations. 

We mention that the scale independence of the un- 
renormalised Yukawa coupling ha seen in figure 10 is also 
found in the present truncation for Up = 0. For other val- 
ues of the parameters the unrenormalised Yukawa cou- 
pling indeed changes more substantially. In the present 
truncation the fiow of TOq is infiuenced only through ha 
(see (|B12|) ). Therefore, if ha is altered during the fiow 
this will result in a change of the final value of Wia and 
hence the critical temperature. It is precisely the inter- 
play of the running of ha with other effects that leads 
to the robustness of the flow equation result. In a sense, 
the flow finds the "correct value" of ha independently 
of the initial value for fc — > oo. Indeed, the initial value 
depends strongly on the choice of Um and this "mem- 
ory" must be lost for physical quantities to be indepen- 
dent of the choice of bosonisation. We also observe that 
the "Schwinger-Dyson value" Up — Um is not the choice 
which minimizes the effect of the bosonic fluctuations. 
This rather happens for Up/t « 2.5. 



The parameterisation independent Hartree-Fock results corre- 
spond to a critical temperature of approximately Tc ~ 2.9£ in 
this case. 



XVII. CONCLUSIONS 

The phase diagram of a high temperature supercon- 
ductor shows many complicated features. At low dop- 
ing and low temperature these materials are antiferro- 
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magnetic insulators. Increasing the concentration of elec- 
trons or holes turns them into a superconductor with 
exceptionally high transition temperatures as compared 
to "conventional" superconductors. Not much is known 
about the mechanism for the binding of electrons into 
Cooper pairs in these materials. Between doping concen- 
trations leading to antiferromagnetic or superconducting 
behaviour one suspects a region in which many different 
degrees of freedom seem to play a role. The clarification 
of the basic degrees of freedom and their interplay in this 
pseudogap region still needs a lot of experimental and 
theoretical effort. 

The common feature of all high temperature supercon- 
ductors is their highly anisotropic structure composed 
of layers of copper oxide (CUO2) planes. The interest- 
ing properties of these materials and the mechanisms for 
generating them seem to be largely confined to these 
planes. The two dimensional Hubbard model is a sim- 
ple attempt to capture this microscopic structure. The 
model assumes electrons that are able to tunnel from site 
to site on a lattice and feel a mutual screened Coulomb 
repulsion. Whether such an oversimplified model is able 
to reproduce the complex phase structure of a real high 
temperature superconductor or parts of it still has to be 
clarified. A lot of theoretical work has been dedicated to 
this task over the last years but so far the results are still 
inconclusive. 

The aim of this paper is a computation of the phase 
diagram, the temperature dependence of the antiferro- 
magnetic order parameter and the critical behaviour of 
the correlation length for the two dimensional Hubbard 
model with next neighbour coupling and small doping. 
Although the quantitative results of this mainly theoret- 
ical investigation cannot be directly compared with ex- 
perimental results we believe that several universal fea- 
tures actually hold independently of the details of the 
superconducting material. We attempt a solution of the 
Hubbard model by means of renormalisation group (RG) 
equations. Earlier RG studies have already revealed the 
power of this technique in the context of the Hubbard 
model 0-01. The use of a purely fermionic language, 
however, renders an investigation of the low temperature 
phase extremely difficult. We believe that it is favourable 
to include the interesting degrees of freedom more explic- 
itly. This can be achieved by rewriting the original action 
of the Hubbard model in a partially bosonised form. 

Already a simple mean field calculation in the partially 
bosonised Hubbard model leads to very encouraging re- 
sults. We were able to reproduce Q a phase diagram that 
closely resembles the one of a real high temperature su- 
perconductor. However, this simple approach also reveals 
an undesirable drawback of the bosonisation procedure. 
The couplings are not uniquely fixed by the reformula- 
tion procedure but there is an arbitrariness connected 
to different parameterisations of the coupling term. This 
"Fierz ambiguity" is also mirrored in the results. Even 
though the reformulation itself is exact, approximations 
break this parameterisation invariance and often lead to 



a strong dependence of the computed phase diagram on 
unphysical parameters. 

In the mean field approximation the fluctuations of 
the bosonic fields are completely neglected. Taking the 
bosonic fields into account should dispose of or at least 
diminish the parameterisation dependence of the results. 
An inclusion of the effective bosonic fluctuations may 
be performed using renormalisation group equations. We 
base our investigation on truncations of an exact func- 
tional renormalisation group equation for the effective 
action. The partially bosonised theory then serves as a 
starting point for the flow of couplings. A one loop calcu- 
lation serves as a guide for the formulation of suitable 
truncations schemes and also clarifies the relation be- 
tween diagrams in the bosonised theory and the original 
fermionic formulation. 

A first truncation deals with antiferromagnetic be- 
haviour at and close to half filling. We are able to ob- 
serve the breaking of the spin rotation symmetry and 
may follow the flow further into the ordered phase. We 
obtain a plausible explanation of why antiferromagnetic 
behaviour may be observed in the two dimensional model 
despite of the Mermin- Wagner theorem which forbids the 
breaking of a continuous symmetry in two dimensions for 
T > 0. The antiferromagnetic order is lost only when av- 
eraging over extremely large spatial extensions, much be- 
yond any practical macroscopic scale. (In a strict sense, 
observation of antiferromagnetism may be regarded as a 
finite size effect.) We emphasize that effective long range 
order is a two-dimensional phenomenon which does not 
need the existence of a third dimension in realistic ma- 
terials (with weak coupling between the plans). For a 
low doping concentration we calculate a phase diagram 
that agrees well with other investigations. 

A necessary check of the consistency of our truncation 
scheme concerns the question to which extend physical 
observables depend on the unphysical parameters intro- 
duced in the bosonisation procedure. We have studied the 
dependence of the pseudocritical temperature Tpc on the 
unphysical ratio of Yukawa couplings Um/Up (at fixed 
U — 3Um — Up). Fig. 14 shows that the Fierz ambiguity 
of mean field theory is substantially reduced by the in- 
clusion of the bosonic fluctuations, even in a very crude 
truncation. This result is encouraging and tests along this 
line should be performed for future more elaborate trun- 
cations. 

Our main result is the temperature dependence of the 
antiferromagnetic order parameter Ka (T) shown in fig. ^ 
This determines the gap for the electrons according to 
Aa — ^/2h^K^t. We also have given a detailed descrip- 
tion of the universal critical behaviour for temperatures 
near the critical temperature Tc, including the tempera- 
ture dependence of the correlation length for T > and 



Of course, the third dimension may nevertheless, play a role for 
some realistic cuprates. 
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an estimate of tlie range of validity of the universal be- 
haviour (T < O.STpc). In agreement with other methods 
[sills^ls^ we find that Tc is substantially lower than the 
Hartree-Fock result. 

All our findings point to a simple Landau theory for 
the antiferromagnetic spin waves in the low temperature 
phase 

1 1 _ / 2T \ ^ 

Sl = -Za{T)d,ad.a+ -Xa{T) i^a' - —^naiT)j . 

(136) 

By variation with respect to a{x) one finds the field equa- 
tions for static spin wave configurations in thermal equi- 
librium 

Aa = ^(d^- -^Ka] a. (137) 

The effective parameters Za,Xa and Ha depend on T, 
the chemical potential /x (doping) and the characteristic 
length scale of the configuration or experimental probe 
k~^. They can be directly extracted from the solution of 
the flow equation. Since Za and Aq depend continuously 
on fj. and are positive for fj, = they will remain positive 
in a whole range of fi. We see no sign of instability to- 
wards a spiral state for small nonzero doping 13, '34^, '3^ . 

In view of the self-consistency of our picture concern- 
ing the Mermin- Wagner theorem and the Fierz ambigu- 
ity we believe that our exploration of the low temper- 
ature phase is reliable. The quantitative accuracy can 
be improved by future more extended truncations. This 



should confirm if all essential physical ingredients are in- 
deed taken into account in the present truncation. It will 
also permit a more detailed study of the universal critical 
behaviour near the second order phase transition and an 
extension of the phase diagram to larger doping. 

Already at the present stage our investigation goes 
substantially beyond the results of a Hartree-Fock or 
Schwinger-Dyson calculation. The additional effects of 
the bosonic fluctuations that we have included in the 
present work lowers the critical temperature by roughly 
factor two. This demonstrates the importance of these 
fluctuations even for the rough features of the phase dia- 
gram. A robust picture of the low temperature phase and 
the phase transition can only be reached if the fluctua- 
tion effects of these collective bosonic degrees of freedom 
are taken into account properly. 

Indeed, the functional renormalisation group in the 
partially bosonised formulation presented here has inter- 
esting prospects for understanding further features of the 
Hubbard model. It allows us to investigate directly the 
low temperature phases with nonvanishing order parame- 
ters. Further channels of possible spontaneous symmetry 
breaking can be described by additional bosonic degrees 
of freedom. This concerns, in particular, the boson re- 
sponsible for d-wave superconductivity for large doping 
[3 which is briefly discussed in appendix D. Adding fur- 
ther bosonic channels also improves the resolution of the 
momentum dependence of the four fermion vertex. It re- 
mains to be seen how these methods can contribute to the 
fascinating theoretical challenge of establishing a quanti- 
tative phase diagram of the Hubbard model. 



I 

APPENDIX A: FOUR-FERMION COUPLINGS 



In this appendix we compute the box diagrams shown in fig. 15. 




FIG. 15: Box diagrams generating four fermion interactions. 



They generate four fermion interactions by the renormalisation How, which will later be absorbed into effective bosonic 
interactions by rebosonisation (sect. |X|) . 
We choose the external momenta as 



A^Q, B^Q-K, C^Q\ D^Q' + K 



(Al) 
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such that for the rh boson 

Arbox= ^QQ-{K)vaQ{K)rhQ,{^K) + --- , 

KQQ' (A2) 

rhQ{K)^4,*{Q)SHQ + K). 
For the p-boson the Pauh-matrices t? are omitted. We evaluate \qq' (K) at 

^ ^/ „ , I = n for the a boson 

Q^Q' ^0, and <^ (A3) 
\K — for the p boson. 

As an example, if only the a boson is present, we obtain 

A A, = -t,Y,p-\L)p-\U - L) Pp\-L)[Pp\U - L) - Pp\L)]. (^4) 

L 

For the evaluation of this expression leading to eqs. H12()(l . (|B10ll . (|B17p . IjBlSjl one employs the identities 



(A5) 



APPENDIX B: EXPLICIT FORM OF THE FLOW EQUATIONS 
1. Antiferromagnetism at half filling 

We display here the explicit flow equations for the mass term, order parameter and quartic scalar coupling for fi ~ 0. 
They are obtained by taking appropriate derivatives of eq. Ulllf) . With the definitions ija = —df lnZa (anomalous 
dimension) and '^qfiQ) = X]q ©(^^ — Q^)f{Q) we get in the symmetric (SYM) regime 

^(^^^a- -2^2^ (fc/t)2(l+™2)2 

+ (d.J.)/^a2fcj.3y_^(2,)4 y(o) +eosh^(0) 
- (2 - f]a)ml, 



2 - 77a(l - QVfc^) 

4 r d^q ftanhy(O) , 2y(0) tanhy(O) - 1 ) (B2) 



'8fcr|7_, (2^)H y'(o) 2/2(0) cosh2y(o) 

-2(l-77a)Aa. 

In the broken (SSB) regime we find 

_ 1 ^2-ry,(l-QVfc2) 



2T^ {k/ty \{l + 2\aKaY 



hi r d^q [ tanhy(0 1 ) (B3) 

^ " \a 2fcT3 (27r)2 \ yi^,) ^ cosh' y{Ka) 



" r 4^ (W l(l + 2A...)3+2j 



.,4 r cg2g f tB.nhy{Ka) 2y(^a) tanh y(^,) - 1 ) (B4) 

^"^'^ «8fcr| 1,(2^)2 \ y3(^^) + y2(^JcOsh2j^(«:„) 

2(l-r;a)Aa, 
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where 

y{^a) = ^ ^e^{q) + 2t^hlKa. (B5) 

If k obeys both k < 2ttT and k < it we are able to evaluate the ^^^^ ^^'^ ^ %(! ^ Q^/^^)) = 

fc2T(4-77,)/(8^). 



2. Antiferromagnetism for nonzero doping 

For /i 7^ the flow equations for the couplings and in the bosonic potential can be obtained by appropriate 
derivatives of p32|l with respect to a 

- (2 - ?7a)TO^, 

A2;^2-r;„(l-QVfc') 



(fc/T)2(l+^2)3 



^/5T^U4 r f tanhy(M) , 2y(A.)tanh^(/i)-^(-A.)/^(0) ) (B7) 



'8fcT5 7_^ (27r)2 \ y3(o) r{0)cosh'y{^l) J 

-2(l-7?a)Aa, 



with 



y(M)-(e(g)-A')/(2rfe). (B8) 
For the flow equation of the Yukawa coupling we evaluate eq. (|12U|I with pf{Q) = {iojp + e(q) — fJ-)/t- One finds 

kdkhl^hl7ja+f3l;i^ +f3l['\ (B9) 



/3g) = - 2ht^kd, J2 {PFliQ)PaiQ)} 

Q 



= - ^rnlhl^kd, J2 {PF^iQ)PaHQ)PaHQ + n)} (BIO) 



'k^T 

Q 



with 



(Q + n)2/P+m2 

^2 (o;^ + + ^)2)2(^| + „ ^)2)2 • 



(Bll) 
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3. Antiferromagnetism and charge density waves 

The running of the anomalous dimension rja can be inferred directly from (|117|l . The flow of the masses is 
governed by a fermionic loop (the bosonic loop does not contribute as we have neglected the quartic bosonic coupling) 
and reads 



(Pq ftanhy 1 



y 



cosh^ y J ' 



(B12) 



^ _2 (Firj^yr^T f d^q J l-ytanhy 



(B13) 



where y = 

The running of the Yukawa couplings is governed by 



(B14) 



with the beta-functions 



(B15) 



4'^) = h,kdkJ2PF\L){3h,p-\L) + h^P-\L)} 



(B16) 



and 



4? ^^^^'^ E {^htp-'imr^n L) P^\~L)[Pp\li -L)- Pp\L)] 



+ hlhl{p-\L)p;\n -L) + p-\L)p-\n - l)) x p^\-L)[Pp\ii + p^^l)]}, 



(B17) 



/3f ) ^^kd,J2{{3h:Pa-'{L)+KP^'iL)) P^\-L)[P^\-L) + P^\L)]}. 



(B18) 



In a more explicit form they read 



I"' (fq ~ 
dkK =haT I -—TTjdk 



1 



h 



ak 



{TiTkf + el\\Pak{q 

1 r h 



Ka\ _ ( 

PaA^ V 



Ppk{q) PpaJ 

A 



hpA\ 1 tanhy ( _ hpA_ 



(7rTfc)2 + el { Pakiq)PakiTT - q) Pa^PaK j 4r|y P2 



^T^y yPaK PpA 
—4 

tanhy /i^^ 



(B19) 



with y — Similarly we find for hp 



dkhp =hpT [ ,f l^ dk 



. (27r)2 



{TrTky 



7-2 



71^ 



—2 —2 —2 



[e2 +(^rfc)2]2 ^VPpfe(q) PpA/ \Pak{q) Pa^n cosh' y I PpA PaA 



pk 



<r2 r 



ak 



aA ' 
2 



'pfe 



—4 -r4 

1 ^tanhy 1 ^ J ^ /i^j. ftpA 



^8?^.'^ y co.h' y'VP!^ ' Pf 



pA 



(B20) 



APPENDIX C: FERMIONIC REGULATOR 



of the fermionic part of the flow equation for the bosonic 



Let us take a look at how the fermionic regulator func- 
tion works. For this purpose we investigate the integrand 
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FIG. 16: Fermionic regulator. We plot the function F{qi,q2 ) defined in lEU for T^/t = 5 (left) and Tk/t = 1/5 (right). 



mass (IBIOI). i.e. the function 



F{qi,q2) 



tanhy 

y 



1 



cosh y 



2t{cosqi + cosg2) 



(CI) 

If the flow is properly implemented F{qi,q2) should be 
dominated by momenta close to the Fermi surface for 
small Tk/t, since we only change the size of the shell of 
fluctuations around the Fermi-surface as k is lowered. In 
fig. 16 we plot F{qi,q2) for different values of Tk/t. The 
left part of fig. 16 shows F{qi,q2) for fairly large values 
of Tk while the right part is for low Tk ■ Observe that the 
contribution to the integral comes from narrower regions 
around the Fermi-surface the smaller Tk becomes. This 
was exactly what was intended by the regulator. 



and similar for x y. 

Let us now introduce auxiliary fields B — (p, to, s, c, Cy) 
and add a term quadratic in these fields to the action such 
that the four fermion interaction is just cancelled: 

s[4',r,B] = s[i',r] + ^s,r.S,r,B], (d4) 

AS"™* = ^ < \Up{pi ~ pif + \Um{rrii ~ rfiif 
X ^ 

+Us{s* - s*){s, - Si) + U^[{c*^, - cli){c^i - c^i) 

-^(Pj - Pi){Pi+e^ - Pt+eJ (D5) 

+ i(mj - TOi)(mj+e^ - rrit+ej] + Uy[x -> y] 



APPENDIX D: COOPER PAIRS IN THE d-WAVE 
CHANNEL 

The Hubbard interaction may be decomposed into 
fermion bilinears in many different ways. To display a 
few, we define fermion bilinears corresponding to charge 
density, magnetisation and Cooper pairs in different 
channels 

pix) = = ^rv-^, (Di) 

rn{X) = rhi ^ 

s{X) = = -0,e-0j, s*(X) = s* --0*e-0*, 

where e is the two dimensional completely antisymmetric 
and ix is the unit vector in x-direction. 



tensor (e 



We also define Cy similar to Cx ■ With these definitions we 



may rewrite the interaction term as follows 
and further note the identity 



p1 



-iTO^ 
3 » 



~ 2 
-TO 



3,i 



2 *i 



-PiPi+e^ + TOjTOi+e_^ -|- 2c* ^4 



(D2) 



(D3) 



Restricting the couplings to the range 

Ur > 0, 

Up,U^ >Ux + Uy, 



(D6) 



ensures that the auxiliary fields are Gaussian and can be 
integrated out after a shift of variables. Imposing 



3C/™ -~Up-2Us^U 



(D7) 



furthermore ensures that the resulting four fermion inter- 
action in ASint cancels precisely the one in the original 
action. These conditions thus ensure that the partially 
bosonised partition function 



Z[n:V*]= j V{TP*,Tp,B)exp(^-S[Tp,'ip*,B]+'n*Tp+'qTp 

(D8) 

is indeed equivalent to 

Collecting the terms in ljD4|) we see that we are now 
dealing with a theory of fermions coupled to bosons 
via a Yukawa interaction. In Fourier space the partially 
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bosonised action reads 

S[i^,r ,B]^ Sk^n['4^,r ,B] + SY[i>A\B], (D9) 

Sktn = X! {'^*('3)Nq - /-^ + 
Q 

+\{Up ~ Uxcosqx - Uy cos qy)p{-Q)p{Q) 
+^{Um + UxCOSQx + Uy COS qy)m{-Q)m{Q) 
+UsS*{Q)s{Q) + U^cl{Q)c^{Q) + C/j,c;(Q)cj,(Q)} 

where for next neighbor [t) and diagonal neighbor (t') 
couphngs one has 

e{q) = —2t{cosqx +cosqy) — it' cos q^ cos qy]. (DIG) 

The Yukawa interaction reads 



Sr = - E 

KQQ' 



5{K-Q + Q') 

{Up - U^cosk^ - Uy cos ky)p{K)^*{Q)TP{Q') 

+ {Ur,^ + COS + Uy COS ky)rh{K)i}* {Q)ai>{Q') 
+5{K -Q-Q') 
Us[s*{K)i^{Q)e^{Q') - s{K)r{Q)er{Q')] 



+ U-r COS 



2 



(i^)^(Q)e,A(Q') 



-c,{K)r{Qyr{Q')] 

-Uy COS ^ [c;(ir)i/;(Q)e^(g') 

-cy{K)r{Q)4*{Q')] 



(Dll) 



Here we have used the Fourier transforms l@J for the 
fermions and 



r(Q) 



for x'^*^ = (p, TO, s'*^), while for c, c* 



we use: 



Q 



(D12) 



(D13) 



(*) 
V ■ 



and similar for c. 

In order to get an operator that has d-wave symmetry 
we perform the transformation of variables 



e = {cx + c.y), c^ = i(e + d), 



(D14) 



and similar for e*,d* . If we insert this variable transfor- 
mation into the action (|D9|I and put Ux = Uy = Uc we 



A 



FIG. 17: The (7 = mode of d(Q) in its local form. Solid and 
dashed lines indicate that two fermion operators on neigh- 
bouring lattice sites are connected with positive or negative 
sign respectively. 



obtain for the e and d dependent part of the action: 

ge,d ^ \U,{e*{Q)e{Q) + d*{Q)d{Q)} (D15) 



5{K - Q - Q')[^(cos 



COS 



KQQ' 



x[e*(if)^(Q)eV'(Q') - e{K)r{Q)4*{Q')] 



+ %(cos- 



— COS ■ 







><[d*{K)ij{Q)ti^{Q')-d{K)r{Q)tr{Q')]]- 

Integrating out the bosons is equivalent to inserting the 
solutions of their field equations 6S[iIj,i/j* , B]/SB — 
(i.e. the saddle point) into the action Slip, iIj*,B]. For the 
boson d the solution is 



diQ) = c.AQ)-CyiQ)^d{Q), 



(D16) 



as expected from the construction of d. In the fermionic 
theory it is thus this combination that is represented by 
the boson d. 

Let us take a look at the q — (0, 0) mode of d{Q), i.e. at 
a spatially homogeneous field. From (|D16p we know that 
it is a superposition of stripes along the x- and y-axis 
added with opposite signs. To find a "local" expression 
we rewrite 



d{Q = Q) = Y{coskx - cosky)tP{K)e4!{-K) (D17) 

K 

-^iX)e^PiX + iy) ~ ^P{X)e^iX - Cy)} 



X 



so that at each lattice site we find an operator of the 
form shown in figure 17. We see that indeed this boson 
may serve as a lattice representation of a d^^^-y^ wave as 
it changes its sign under rotation by 90° but not under 
reflection at the x 01 y axes (see also [s^l for a more ex- 
tensive classification). The mean field approach including 
the boson d is discussed in 
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